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ABSTRACT 


At  Reynolds  numbers  greater  than  about  5  x  10  ,  corresponding  to 
altitudes  below  about  l8o, 000  ft.,  the  hot  "outer  inviscid  wake" 
behind  the  bow  shock  wave  produced  by  a  blunt-nosed  body  at  hypersonic 
speeds  is  cooled  mainly  by  turbulent  diffusion  and  conduction.  Turbu¬ 
lence  originates  in  the  "inner  wake"  formed  by  the  coalescence  of  the 
free  shear  layers  (or  annulus)  shed  from  the  body  surface  when  the 
boundary  layer  separates  from  the  surface.  As  this  turbulence 
spreads  outward,  it  swallows  enthalpy  or  momentum  defect  originally 
contained  in  the  outer  inviscid  wake.  If  the  turbulence  is  "locally 
similar",  i.e.,  if  it  behaves  at  each  station  like  a  slice  of  a  low 
speed  "self -similar"  wake,  then  the  turbulent  diffusivity  grows  from 
a  low  initial  value  near  the  body  to  a  value  corresponding  to  the 
total  drag  of  the  body  at  about  300  body  diameters  downstream.  At 
flight  velocities  of  the  order  of  9, 000-10, 000  ft. /sec.  the  growth  of 
the  turbulent  inner  wake  predicted  on  the  basis  of  "locally  similar" 
turbulence  is  in  good  agreement  with  shadowgraph  measurements  of 
wake  widths  behind  spheres  obtained  in  ballistic  ranges  in  the 
region  from  200  to  4,000  body  diameters  downstream  of  the  body. 
Tentatively,  one  concludes  that  the  turbulence  mechanism  in  the  wake 
with  respect  to  a  fixed  observer  is  similar  to  the  low  speed  case,  in 
spite  of  the  large  mean  temperature  gradients. 


electron  density  in  a  turbulent  wake  "behind  a  "blunt,  "body,  the  two 
limiting  cases  of  thermodynamic  equilibrium  and  pure  diffusion  (zero 
electron-ion  recombination  rate)  are  calculated  for  =  22  at 
altitudes  of  100,000  ft.  and  200,000  ft.  Even  for  the  case  of  thermo¬ 
dynamic  equilibrium,  the  predicted  turbulent  radar  trail  length  is 
about  200  body  diameters  at  L-band  (1300  mc/sec.)  at  100,000  ft. 
altitude  and  about  15 0  body  diameters  for  UHF  (400  mc/sec.)  at 
200,000  ft.  One  interesting  result  is  that  the  width  of  the  plasma 
cylinder  corresponding  to  the  plasma  frequency  at  L-band  remains 
virtually  constant  at  about  3*5  body  diameters  in  the  range  30  x/d  <. 
150  at  100,000  ft.  altitude.  These  results  are  sufficiently  encourag¬ 
ing  so  that  one  can  think  about  including  the  effects  of  finite 
chemical  and  electron  -  ion  recombination  rates  in  the  analysis  in 
order  to  give  a  more  complete  picture  of  the  wake  at  hypersonic  speeds. 
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1.  Introduction 

The  wake  produced  by  a  body  moving  through  a  fluid  medium  is  one  of 
the  oldest  known  fluid-mechanical  phenomena.  In  calm  water  the  wake  of 
a  sailing  ship  is  observed  for  many  miles  behind  the  ship.  Experi¬ 
mental  studies  of  low  speed  wTakes  show  that  the  momentum  defect 
contained  in  the  viscous  portion  of  the  wake  spreads  rather  slowly  into 
the  rest  of  the  fluid,  even  when  the  diffusion  process  is  turbulent, 

High  speed  wake  phenomena  have  also  been  noticed  for  a  long  time,  but 
much  less  is  known  about  them.  When  the  familiar  meteor  trails  are 
studied  by  radar  echo  techniques,  trails  at  least  15  miles  long  are 
recorded  at  altitudes  of  the  order  of  100  miles.  For  small  meteors 
the  trail  length  is  determined  partly  by  the  time  required  to  consume 
the  object,  and  partly  by  laminar  diffusion  and  electron  attachment 
processes  in  the  wake. 

Wakes  produced  by  slender  projectiles  and  blunt  bodies  at  supersonic 
speeds  have  been  observed  in  ballistic  ranges  by  shadowgraph  and  other 
optical  techniques,  but  until  recently  these  observations  did  not  extend 
very  far  behind  the  body.  One  of  the  most  striking  features  of  the  flow 
field  around  a  blunt-nosed  object  at  supersonic  speeds  is  the  bow  shock 
wave  (Figure  l).  When  a  missile  or  spacecraft  enters  the  earth’s 
atmosphere  at  hypersonic  speeds  a  distinct,  nearly  paraboloidal  shock 
surface1*  is  first  formed  in  front  of  the  body  at  altitudes  of  the 
order  of  350,000  -  400,000  feet.  The  gas  that  traverses  the  "strong", 
nearly-noxmal  portion  of  the  bow  shock  wave  is  compressed  and  heated 


*  Superscripts  denote  references  listed  at  end  of  paper. 
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irreversibly  in  the  shock,  and  forms  an  "outer  wake"  behind  the  body 
(Figure  1).  Even  after  isentropic  expansion  back  to  ambient  pressure, 
which  occurs  in  a  distance  of  50  -  100  diameters  behind  the  body,  the 
maximum  gas  temperature  in  this  outer  wake  in  thermodynamic  equi¬ 
librium  is  about  4500°K  at  satellite  velocity,  and  about  6500°K  at  a 
flight  velocity  of  40,000  ft. /sec.  .  Clearly  the  trail  length  defined 
in  terms  of  observables  such  as  electron  density  or  gas  radiation 
intensity  is  determined  by  the  diffusion  and  recombination  processes 
in  this  cylinder  of  hot  gas  left  behind  in  the  atmosphere. 

2 

Feldman  treated  this  problem  under  the  assumption  of  thermodynamic 
equilibrium  and  a  laminar  or  molecular  diffusion  process.  Under  these 
assumptions  the  problem  downstream  of  the  "pres sure -controlled"  portion 
(x/d/>50)  is  closely  analogous  to  the  cooling  of  a  gas  cylinder  with 
a  given  initial  temperature  distribution.  The  time  required  to  reduce 
the  enthalpy  is  proportional  to  r^  ,  where  K.  =  k/^>  cp  is  the  thermo¬ 
metric  conductivity  and  rn  is  the  body  nose  radius,  which  characterizes 
the  breadth  of  the  initial  enthalpy  distribution.  Here  time  is  equiva¬ 
lent  to£/u  ,  wherein  is  distance  behind  the  body,  so  that 
£/r~u  r  /r  =  Be  x  Pr  where  He  is  the  Reynolds  number  and  Pr  the 
Prandtl  number.  Thus  the  predicted  laminar  trail  length  is  very  long 
even  at  high  altitudes;  for  example  for  a  nose  radius  of  one  foot  at 
250, 000  feet  altitude  £/rn  ^  10  -  10  ,  depending  on  the  observable  of 

interest  while  at  60,000  feet  altitude  A/r  '--'10'’  -  106. 

’  '  n 


V 
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Feldman  himself  suggested  that  the  cooling  process  io  much  faster  if 
turbulence  is  produced  in  the  vake .  One  expects  turbulence  to  originate 
in  the  flow  region  with  the  highest  velocity  gradients,  and  this  region 
is  undoubtedly  the  "inner  vake"  formed  by  the  coalescence  of  the  free 
shear  layers  (or  annulus)  shed  from  the  body  surface  when  the  boundary 
layer  leaves  the  surface  (Figure  l).  Experimental  evidence^  shows  that 
transition  to  turbulent  flow  occurs  downstream  of  the  "neck"  at  fairly 
low  Reynolds  numbers,  corresponding  to  altitudes  of  the  order  of  180,000 
feet.  Near  the  neck  this  turbulence  is  confined  to  a  narrow  region 
around  the  wake  axis,  but  it  soon  spreads  outward  by  feeding  on  the 
surrounding  gas.  Eventually  all  of  the  streamlines  originally  in  the 
"outer  wake"  are  engulfed  by  the  turbulent  inner  wake.  The  main 
purpose  of  this  paper  is  to  analyze  this  swallowing  process,  which 
determines  the  rate  of  diffusion  and  cooling  in  the  wake  of  a  blunt- 
nosed  body  at  hypersonic  speeds. 


In  order  to  bring  out  the  main  ideas  as  simply  as  possible  the 
assumption  of  theimodynamic  equilibrium  is  retained  in  the  present 
paper.  The  influence  of  finite  recombination  rates  can  be  included 
later  when  theVturbulent  diffusion  process  is  somewhat  clearer. 


In  Section  2  the  structure  of  the  wake  behind  a  blunt-nosed  body 
at  hypersonic  speeds  is  described,  and  simplified  representations  of 
the  outer  and  Inner  wakes  are  introduced.  The  boundary  between  these 
two  regions  is  supposed  to  be  a  sharp  '''front"’,  and  the  growth  of  the\^  j 

y 


-t<>  /> 


inner  wake  depends  only  on  the  gradient  and  value  of  the  enthalpy  at 

this  front.  Two  limiting  cases  of  the  behavior  of  the  turbulent 

diffusivity  are  studied:  (l)  ^locally  similar”  turbulence,  in  which 

the  flow  at  each  station  behaves  like  a  "slice”  of  a  low  speed  "self- 

similar"  turbulent  wake,  so  that  the  diffusivity  is  proportional  to 

the  local  momentum  defect  or  dmg  contained  in  the  inner  wake;  (2) 

T  i 

frozen  diffusivity,  in  which  the  turbulent  diffusivity  depends  only 
on  the  initial  value  of  the  drag  coefficient  for  the  inner  wake  at  the 
' neck  . v  Once  the  diffusivity  is  specified  the  turbulent  diffusion 
equations  for^enthalpy  and  mass  concentration  can  be  integrated 
(Section  3)*  In  Section  4, ^specific  examples  of  the  growth  of  the 
turbulent  inner  wake  are  calculated  for  comparison  with  measurements 
made  in  ballistic  ranges.  A  typical  reentry  example  is  also  computed 
at  Mcg^=  22  and  an  altitude  of  100,000  feet,  in  order  to  illustrate  the 
behavior  of  the  enthalpy  distribution,  and  the  behavior  of  the  profiles 
of  electron  density  in  the  wake  in  the  two  limiting  cases  of  thermo¬ 
dynamic  equilibrium  and  pure  diffusion  (zero  recombination). ^Conclusions 
about  the  hypersonic  wake  obtained  from  this  study  are  summarize aln. 
Section  5,  which  also  contains  some  suggestions  for  future  work. 
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2.  Model  of  the  Hypersonic  Wake 

2.1.  Description  of  the  Flow  Pattern 

In  order  to  determine  the  c 02*1*6 ct  initial  conditions  for  the 
lateral  spreading  of  the  turbulent  inner  wake,  we  must  examine  the 
main  aspects  of  wake  formation  behind  a  blunt  body  at  hypersonic  speeds. 
Shadowgraph  observations  supplemented  by  simple  physical  considerations 
show  that  the  hypersonic  viscous  wake  differs  from  its  low  speed 
counterpart  in  three  important  respects:  (l)  no  evidence  of  large-scale 
vortex  formation,  even  for  laminar  flow;  (2)  much  greater  "stability” 
of  the  laminar  free  shear  layers;  (3)  an  initial  value  of  the  momentum 
defect  or  "drag"  contained  in  the  inner  wake  (Figure  1)  one  to  two  orders 
of  magnitude  smaller  than  the  total  body  drag. 

Even  on  a  bluff  body,  such  as  a  sphere,  separation  occurs  well  aft 
of  the  90°  station,  partly  because  of  the  falling  pressure  over  the 
front  half  of  the  body,  and  partly  because  the  supersonic  shear  layers 
can  negotiate  relatively  high  reccmpression  pressure  ratios  (2-3)  at 
the  "neck"  (Figure  l).  The  free  shear  layers  (or  annulus)  shed  from 
the  separation  line  (or  ring)  on  the  body  surface  show  a  "rigidity" 
characteristic  of  supersonic,  and  especially  hypersonic  flows.  Large- 
scale  vortex  formation  just  behind  the  body  is  apparently  inhibited  by 
the  fact  that  violent  streamline  deflections  would  be  accompanied  by 
large  pressure  changes  and  shock  wave  fomation.  The  separated  flow 
tends  to  retain  the  direction  established  at  the  point  of  separation. 


and  the  flow  converges  toward  a  minimum  section  or  "neck",  which 
contains  the  effective  rear  stagnation  point. 

When  the  boundary  layer  on  the  body  surface  is  turbulent  of  course 

the  free  shear  layer  (or  annulus)  is  also  turbulent;  however,  if  the 

boundary  layer  is  laminar,  the  laminar  free  shear  layers  show  a 

remarkable  persistence  at  supersonic  speeds.  Experimental  Btudies 

(and  some  stability  considerations)  indicate  that  the  transition 

Reynolds  number  Re^,  for  the  laminar  free  shear  layer,  based  on 

properties  at  the  outer  edge  and  the  length  of  laminar  "run",  increases 

4 

very  rapidly  with  increasing  local  Mach  number  .  For  an  insulated 

^  k  _ _ _  tr 

body  Re  =  5  x  10  at  M  -  0.5,  Re  =  2  x  lCr  at  M  =  2.0,  and 
Re^  =  2-3  x  10^  at  M^s  4.5,  the  highest  Mach  number  studied.  Partly 
this  effect  is  the  result  of  the  high  average  kinematic  viscosity  across 
the  layer^ .*  Downstream  of  the  neck,  on  the  other  hand,  the  relevant 
Mach  number  based  on  the  difference  between  the  flow  velocities  on 
the  axis  and  at  the  edge  of  the  inner  wake  decreases  rapidly  with 
distance  away  from  the  neck,  and  the  flow  apparently  becomes  unstable . 

In  addition,  the  inner  wake  is  probably  subjected  to  disturbances 
originating  at  the  neck.  Recent  experiments  at  GALCIT  indicate  that 
the  local  transition  Reynolds  number  in  the  inner  wake  behind  an 


*  This  fact  was  pointed  out  to  us  by  Dr.  Edward  Zukoski,  on  the 
basis  of  his  experience  with  hot  wakes  behind  flame -holders. 
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insulated  cylindrical  rod  transverse  to  the  airstream  at  =  5*8  is  of 

If 

the  order  of  5  x  10  ,  based  on  local  properties  and  the  distance  between 

3 

the  neck  and  the  transition  zone .  These  considerations  suggest  that  the 
regime  in  which  the  free  shear  layer  is  still  laminar  upstream  of  the 
neck,  but  transition  to  turbulent  flow  occurs  not  far  downstream,  is 
important  in  understanding  diffusion  processes  in  the  wake  behind  a 
blunt  body  at  hypersonic  speeds. 

If 

As  explained  by  Chapman  the  fluid  below  the  zero  streamline  origi¬ 
nating  at  the  separation  point  on  the  body  is  turned  back  at  the  neck 
because  of  the  pressure  rise  produced  by  the  flow  deflection.  This 
fluid  enters  the  zone  of  recirculation  just  behind  the  body.  The  fluid 
along  the  zero  streamline  is  brought  to  rest,  while  the  fluid  above  the 
zero  streamline,  although  slowed  down,  flows  on  to  form  the  inner  wake 
(Figure  1).  In  Chapman's  limiting  case  the  velocity  along  the  zero 
streamline  just  upstream  of  the  neck  is  about  60  per  cent  of  the 
velocity  at  the  outer  edge  of  the  free  shear  layer.  Thus  the  static 
enthalpy  at  the  outer  edge  of  the  viscous  wake  is  about  35-40  per  cent 
of  the  stagnation  enthalpy  in  the  inviscid  flow.  This  is  to  be  compared 
with  the  enthalpy  along  the  zero  streamline  from  the  free  shear  layer 
which  if  brought  to  rest  adiabatically  has  an  enthalpy  of  about  60  per 
cent  of  the  stagnation  enthalpy  in  the  inviscid  flow,  when  the  body 
surface  is  "cold".*  Thus  a  hot  core  of  fluid  is  generated  in  the  inner 


*  If  the  body  surface  is  insulated,  the  enthalpy  on  the  zero  stream¬ 
line  is  about  95  per  cent  of  the  total  enthalpy  outside  the  viscous  wake. 
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wake,  but  this  core  cools  off  with  distance  downstream  of  the  neck. 
Typical  velocity  and  enthalpy  profiles  just  upstream  and  downstream  of 
the  neck  are  shown  schematically  in  Figure  2. 

The  initial  momentum  defect  or  drag  contained  in  the  inner  wake 
represents  the  sum  of  the  skin-friction  drag  on  the  body  and  the 
momentum  defect  associated  with  the  pressure  rise  at  the  neck.  Since 
the  density  of  the  fluid  in  the  inner  wake  just  downstream  of  the  neck 
is  much  lower  than  ambient  and  the  neck  width  is  usually  less  than  one- 
half  the  body  diameter,  this  initial  drag  is  only  a  small  fraction  of 
the  drag  associated  with  the  bow  shock  at  hypersonic  speeds .  As  we 
shall  see  later  (Section  2.3)  the  turbulent  eddy  diffusivity  £T  is 
proportional  to  the  drag  coefficient  of  the  inner  wake .  We  conclude 
that  the  initial  value  of  CT  is  much  lower  than  the  final  value 
(x/d — corresponding  to  the  total  body  drag.  This  fact  has 
important  consequences  for  the  whole  process  of  turbulent  diffusion 
(Sections  2.3  and  3). 

Near  the  neck  the  amount  of  momentum  defect  swallowed  up  by  the 
growing  inner  wake  is  still  small  compared  with  the  initial  momentum 
defect,  and  the  inner  wake  develops  in  the  downstream  direction  almost 
as  if  it  were  enveloped  by  an  external  environment  uniform  in  the 
direction  normal  to  the  wake  axis.  In  addition  to  the  usual  spreading 
of  the  wake  by  "mixing",  the  physical  boundaries  of  the  wake  expand 
because  the  static  pressure  is  falling  in  the  stream  direction,  and 
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the  whole  flow  field  is  expanding.  The  proper  length  scale  for  the  initial 

development  of  the  inner  wake  is  not  the  body  diameter,  hut  the  momentum 

thickness  0  of  the  inner  wake,  which  is  at  least  ten  times  smaller.  Thus 

the  velocity  difference  Au  =  juf  -  u  (O)]  across  the  inner  wake  decays 

to  a  small  fraction  of  (u»)  .  in  a  downstream  distance  of  the  order 

f  neck: 

of  100  0  to  200  0,  which  corresponds  to  about  10  -  20  body  diameters, 

at  most .  In  fact  when  uf  ~  1  the  mean  inner  wake  flow 

Uf 

downstream  of  transition  to  turbulence  should  exhibit  a  similarity  of 
the  form 

-rtv.4 1 
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_  _ 

— 1 

X  -xc 
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just  as  Townsend  found  for  the  low-speed  wake  behind  a  cylindrical  rod 
transverse  to  the  air  stream.*  (Here  Y  is  essentially  the  stream  function, 
and  m  is  a  geometric  index;  m  =  0  for  a  two-dimensional  flow  and  m  =  1 
for  axially-symmetric  flow.)  A  typical  velocity  distribution  at  x/d"^ io 
is  sketched  in  Figure  3a;  and  the  corresponding  enthalpy  profile  appears 
as  a  "turbulent  pulse"  at  the  base  of  the  "inviscid"  profile  of  the 
outer  wake  (Figure  3^). 


When  x/d  is  of  the  order  of  10,  the  rate  at  which  the  inner  wake  is 
swallowing  momentum  defect  originally  contained  in  the  outer  wake  can  no 


*  Except  possibly  for  a  small  effect  caused  by  the  negative  pressure 
gradient  along  the  wake  axis. 
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longer  "be  ignored.  In  the  limiting  case  of  ''frozen"  diffuslvity  the 

turbulent  eddy  diffusivity  6  ^  retains  the  value  established  by  the 

initial  drag  of  the  inner  wake  at  the  neck.  But  in  the  opposite 

limiting  case  of  "locally  similar"  turbulence,  the  local  value  of 

6^(x)  depends  on  the  local  drag  contained  in  the  inner  wake  at  that 

station.  Thus  (d  £  ^/&x)  depends  on  the  rate  of  increase  of  momentum 

defect  in  the  inner  wake,  which,  in  turn,  depends  on  the  product  of 

the  slope  of  the  inviscid  enthalpy  profile,  (dh/by)  ,  and  the 

y=yf 

rate  of  spreading  (dy^/dx)  of  the  inner  wake.  But  (dy^/dx)  itself 

depends  on(~„.  When  y-/d  <<  1,  (bh/b  y)  is  small,  and  and 
1  x  y-y  -p  ^ 

yf  both  increase  slowly  with  distance  downstream.  Soon  the  turbulent 
front  reaches  the  portion  of  the  inviscid  enthalpy  profile  where 
(bh/by)  is  appreciable,  and  one  expects  a  rapid  growth  ofc^  end 
inner  wake  width.  Eventually  (bh/by)  _  decreases  again  and  the 

y-yf 

rate  of  growth  of  the  Inner  wake  slows  down.  Far  downstream  (x/d  >77  1) 
most  of  the  wake  is  turbulent;  the  turbulent  diffuslvity  takes  on  a 
value  corresponding  to  the  total  drag  of  the  body,  and  similarity  is 
again  established. 


2.2.  Characteristics  of  the  Outer  and  Inner  Wakes 
2.2.1.  Outer  Wake 

In  order  to  bring  out  the  main  features  of  the  px'ohlem  as  simply 
as  possible,  we  consider  Reynolds  numbers  large  enough  so  that  the 
time  scale  for  laminar  diffusion  in  the  outer  wake  is  much  longer  than 
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the  time  required  for  the  turbulent  inner  vake  to  swallow  most  of  the 
momentum  defect.  The  flow  at  any  point  in  the  outer  wake  is  supposed 
to  he  unaffected  by  the  inner  wake  until  the  turbulent  "front"  reaches 
it.  Thus  the  enthalpy  profile  in  the  outer  wake  is  completely  specified 

-ft 

by  the  inviseid  flow  field.  The  present  analysis  is  concerned  with 
the  region  x/d  }  2  -  10,  where  the  static  pressure  is  only  about  5  times 
ambient  pressure,  at  moat.  In  this  region  the  enthalpy  along  an  inviseid 
streamline  is  a  slowly  varying  function  of  static  pressure,  and  we  shall 
ignore  the  streamwise  enthalpy  gradient,  but  evaluate  the  local  enthalpy 
correctly.  This  procedure  amounts  to  a  kind  of  local  similarity  approxi¬ 
mation.^ 


As  in  most  compressible  flows  the  "proper  distance"  normal  to  the 
wake  axis  is  not  the  physical  distance  y,  but  the  Howarth-DorodnitSyn 
variable  y  defined  by  the  relation 

YiT  jy,.  =  (Mv>  .  a) 

where  y  =  y/d,  and  the  subscript  "L"  denotes  "inviseid"  or  outer  wake 
quantities.  The  inviseid  enthalpy  profile  is  represented  by  the  two 
parameter  relation 


(WO _l  -Ovk-J-1.,  3  (°vYu)  = 


(2) 


*  The  wake  shock  (Figure  l)  introduces  a  slight  modification 
(Appendix  l). 
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\ 

for  YT  >  Yr  ,  where  h  denotes  the  complete  static  enthalpy,  including 
-L  Lri 

X  r> 

the  chemical  enthalpy  ,  (h_ )  is  the  value  that  the  enthalpy  would 

Li  o 

have  along  the  zero  streamline  if  the  inviscid  flow  extended  to  the 
hL 

axis,  H  =  (r —  -  1)  ,  z  =  o(  Yt  ,  and  oi  is  a  scale  factor.  In  the 
n  o  m  Li  m 

11  co  \x  ^ 

entire  hypersonic  wake  —  ■  ■  - -  «.  1  and  h  -  h  ^  -u  (u-u  ). 

u  co  ao  <X>  GO 

To  this  approximation  YT  (Eq.  (1)  )  represents  the  stream  function,  or 

X j 

mass  flux  between  the  axis  and  y,  and  the  parameter  cl  is  determined 

m 

by  the  momentum  defect  in  the  inviscid  flow.  In  the  region  x/d}5  -  10 
the  drag  contribution  made  by  the  static  pressure  is  negligible  compared 
to  the  momentum  defect,  and  (Appendix  l), 


Ttt+l 

~rt\ 
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"Do  VU™  '  co 
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(A.l-9) 


Here  is  the  "inviscid"  drag  coefficient  and 
o 


= 


aO 

f  tn. 

i. 


(A.  1-10) 


When  x/d  ^  50  -  100  the  static  pressure  is  virtually  equal  to 
the  ambient  pressure,  the  streamlines  in  the  inviscid  flow  are  all 
parallel  to  the  wake  axis,  and.  the  inviscid  enthalpy  profile  described 
by  Eq.  (2)  is  independent  of  axial  distance,  in  the  limiting  cases  of 
"frozen"  flow  or  thermodynamic  equilibrium.* 


*Evidently  the  methods  of  the  present  paper  can  be  generalized 
to  include  chemical  and  electron-ion  recombination. 
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2.2.2.  Inner  Wake 

The  present  study  is  simplified  "by  considering  only  the  thermo¬ 
dynamic  effect  of  the  static  pressure  variation  along  the  wake  axis, 
hut  not  its  dynamic  effect.  Thus  the  shape  of  the  turbulent  mean 
enthalpy  profile  across  the  inner  wake  is  invariant  in  the  proper 
coordinates  for  a  short  distance  downstream  of  the  neck,  and.  this 
same  "similarity  profile"  should  again  appear  far  downstream.  Although 
the  profile  shape  in  the  intermediate  region  undoubtedly  changes  to 
some  extent,  the  simplest  assumption  is  that  the  mean  enthalpy  is 
described  completely  by  two  parameters  —  the  amplitude  (h^  (0)  -  hf ) 

and  the  wake  width  Y  ,  i.e., 

if 


Since  the  local  environment  surrounding  the  inner  wake  is  characterized 
by  the  values  of  h^,  and  u^  determined  by  the  inviscid  flow  at 
y  =  yf>  the  Howarth-Dorodnitsyn  variable  Y^  for  the  inner  wake  is 
logically  defined  as  follows: 

T-T  JYt  -  %  ^  (4) 

*  A  more  refined  analysis  taking  the  change  in  profile  shape 
into  account  can  always  be  carried  out  once  the  turbulent  diffusion 
process  is  understood  somewhat  better  than  at  present.  Also,  the 
function  G(YT/Yrr  )  is  arbitrary  and  for  the  present  calculations  is 

if 

taken  as  a  parabola. 
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(The  relation  "between  the  inviscid  and  turbulent  mass  flux  quantities 

Yt  and  Y_  is  discussed  in  Appendix  2.) 

Lf  Tf 

x  -  x.  o 

The  two  unknown  functions  ^ (  j  )  and  Y^  (  ;  ) 


are  to  be  determined  by  satisfying  the  following  two  conditions  (Section  3) 
(1)  conservation  of  energy,  as  expressed  by  the  energy  integral  equation 
across  the  inner  wake;  (2)  energy  balance  along  the  inner  wake  axis, 
which  relates  the  rate  of  cooling  in  the  stream  direction  to  the  rate 
of  turbulent  heat  conduction  normal  to  the  axis.  In  the  absence  of  any 

o 

acceptable  theory  of  turbulent  transport  we  adopt  Reynolds'  hypothesis 
of  similarity  between  the  turbulent  transfer  of  mass,  momentum  and 
energy.  Reynolds'  analogy  states  that 


-V  U. 


vVt!.  __  ^ 


(die 


=.  C 


T 


(5) 


where  h^  is  the  thermal  enthalpy  =  I  c^  dT,  is  the  mass  fraction  of 
"fell  ^ 

the  i  species,  and  ^  ^  is  the  common  turbulent  "exchange  coefficient" . 
In  other  words  the  turbulent  Prandtl  and  Lewis  numbers  are  both  equal 

7 

to  unity,  and' 


3 


(6) 


In  his  thorough  experimental  study  of  the  low- speed  turbulent 

£ 

wake  behind  a  long  cylindrical  rod  normal  to  the  airstream,  Townsend 
found  that  a  relation  of  the  form  T  =  p  £.  T  ( c)u/ by)  is  indeed 
satisfied  by  the  mean  velocity  field  not  too  close  to  the  wake  edges, 


2-11 


once  the  mean  flow  becomes  "self -preserving",  or  "similar".*  This  fact 
must  be  regarded  as  an  interesting  coincidence  for  symmetrical  wakes 
(Appendix  3 ) •  As  Townsend  points  out,  pure  gradient  diffusion  cannot 
possibly  describe  the  process  of  envelopment  of  fluid  by  the  large- 
scale  eddies  at  the  outer  boundaries  of  the  wake.  Nevertheless  this 
empirical  description  is  correct  near  the  wake  axis,  and  is  therefore 
extremely  useful  in  satisfying  the  energy  equation  on  the  wake  axis 
(condition  2,  above). 

In  Townsend’s  case  the  flow  outside  the  wake  is  irrotational  and 

the  wake  drag  is  constant.  Once  the  influence  of  the  initial  conditions 

dies  out  the  wake  flow  becomes  "similar",  and  the  wake  is  characterized 

Aul 

by  a  "universal"  Reynolds  number  IL,  =  -  =  constant,  where  Au 

1  €.  T 

is  the  mean  velocity  difference  across  the  wake  and  is  a  measure 
of  the  wake  width,  i.e.,  £  ^  =  (R^)  £o.**  In  spite  of  the 

fact  that  the  inner  wake  exhibits  similarity  at  hypersonic  speeds  only 
in  the  initial  and  final  stages  we  make  the  crucial  assumption  that 
all  along  the  wake  axis 

6.  T  =  K'  Au  yf  ,  (7) 

*  In  Townsend’s  notation  =  \)  ^. 

**  Townsend  found  the  value  Rfp  =  12.5  for  the  wake  behind  a 
cylindrical  rod.  Although  data  on  axially-symmetrie  turbulent  wakes  are 
extremely  scarce  even  at  low  speeds,  the  few  measurements  made  by  Hall^ 
and  Hislop  lead  to  the  provisional  value  R_  ~  14. 
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where  K'  is  a  non-dimens icnal  quantity  dependent  on  fluid  properties  (see 
below),  but  not  on  the  velocity  difference  or  the  wake  width.  According 
to  Eq.  (5)>  this  assumption  is  equivalent  to  the  statement  that  the 
Intensity  of  turbulence  at  each  station  is  proportional  to  th.a  local 
value  of  A  u,  and  the  scale  is  proportional  to  the  local  wake  width 
(Appendix  3)*  Xu  other  words  the  effect  of  previous  history  on  the 
development  of  the  inner  wake  is  ignored,  so  far  as  the  turbulent 
diffusivity  is  concerned.  The  characteristic  time  required  for  internal 
adjustment  of  the  turbulent  motion  is  assumed  to  be  very  short  compared 
to  the  time  scale  for  the  rate  of  change  of  £  in  the  axial  direction, 
and  (=.  T  is  supposed  to  take  on  instantaneously  the  value  appropriate 
to  a  locally  similar  flow.  Later  (Sections  2.3  and  3)  we  shall  also 
consider  the  opposite  limiting  case  of  "slowly  adjusting"  turbulence, 
or  "frozen"  diffusivity,  where  the  quantity  (pf/( lo)  YT  m^T  is  determined 

only  by  the  initial  conditions  in  the  inner  wake. 

The  quantity  K*  could  be  regarded  as  an  empirical  function  to  be 
determined  experimentally,  but  it  is  interesting  to  see  whether  K*  can 
be  related  to  Townsend's  for  low-speed,  isothermal  wakes.  Now  it  is 
well  known  that  the  transformation  given  by  Eq.  (k)  reduces  the  heat 
conduction  problem  with  variable  density  to  an  equivalent  "incompressible" 
heat  conduction  problem.  The  heat  conduction  term  in  the  energy  equation 
specifies  the  equivalent  turbulent  diffusivity.  By  utilizing  Eq.  (k) 
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one  finds  that 

_L  _ _  -^L_  (JLfe  £fT  dk 
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On  the  axis  this  term  "becomes 


But  (Eq.  k), 
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and  the  heat  conduction  term  is  reduced  to  the  form 


<>+') 


VT=0 

where  the  equivalent  turbulent  diffusivity  is  given  "by  the  relation 


*  Following  Mager's**^  work,  Ting^  and  Libby  suggest  a  similar 
relation  of  the  formi^  =  (  p/j5  f)26l  ^  for  two-dimensional  turbulent 

wakes  or  jets  (in  our  notation,  where is  constant).  They  correctly 
point  out  that  such  a  relation  leads  to  a  compressible  eddy  viscosity €.^ 

that  varies  across  the  flow.  In  the  present  study  only  turbulent  heat 
conduction  on  the  wake  axis  is  involved,  so  this  difficulty  does  not 
arise.  Besides,  the  whole  semi-empirical  notion  of  the  existence  of 
6. T  is  probably  correct  only  near  the  wake  axis  (Appendix  3)*  On  the 

wake  axis  Eq.  (8)  also  agrees  with  their  relation  between  ^  and&T 

for  axially-symmetric  wakes  and  jets. 
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According  to  Eqs.  (7)  and  (8), 


then 

6-r  K  All.  Yt  (10b) 


exactly  as  in  Townsend’s  case^j  therefore  K  is  proportional  to 
The  constant  of  proportionality  depends  only  on  the  relationship 
between  our  Y  and  Townsend' s  Jl  (Appendix  3). 

if 


By  adopting  the  value  of  K'  specified  by  Eq.  (10a),  &  T  is  given 
more  precisely  by  the  following  relation  (Eqs.  (7)  and  (8)  ): 


T  _ 


(11) 


When  x/d  J  5  -  10  the  total  enthalpy  is  virtually  a  constant  across  the 
inner  wake,  equal  to  the  total  enthalpy  of  the  inviscid  flow.  Also 
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- -  1  so  that  in  the  turbulent  portion  of  the  flow, 

U<5<5 

=  -ft  -u-co)  (12) 

where  j5uoo  represents  some  average  velocity,  and  O.B  =  |2>  =  0.9«  To 
this  approximation,  Eq.  (11 )  for  £  becomes 


2.3 »  Relation  between  Equivalent  Turbulent  Diffusivity  and 
Momentum  Defect:  Growth  of  Inner  Wake 

2.3.1.  Turbulent  Diffusivity  and  Momentum  Defect 

Once  we  assume  that  the  equivalent  turbulent  diffusivity  is  given 
by  Eq.  (10b)  or  Eq.  (13 ),  it  is  not  surprising  that  the  relation  between 
and  the  local  value  of  the  momentum  defect  in  the  inner  wake  takes 
the  same  form  as  it  does  in  an  incompressible  "similar"  wake  with 
constant  drag.  The  momentum  defect  with  respect  to  the  local  inviscid 
flow  at  the  turbulent  front  is  given  by 


jj3  -  u)  'lj  dij  * 


(Ik) 


*  The  relation  between  this  momentum  defect  and  the  ordinary 
drag  is  discussed  briefly  in  Appendix  4. 


By  utilizing  Eqs.  (3),  (4),  and  (12),  one  finds  that 
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where 
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(15a) 


and 


y  -  v5 


T, 


with  G(  f )  arbitrary  (for  the  present  calculations  in  Sections  2  and 
h  a  parabola  is  chosen).  By  combining  Eqs.  (13)  and  (15)  we  find  that 
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Here  C  '  is  the  "local"  drag  coefficient  defined  by  the  relation 
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where  C  is  the  drag  coefficient  based  on  — - - — 
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.  By  Eq.  (4) 
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(l8a) 


For  "locally  similar"  turbulence  C_  in  Eq.  (l8a)  is  the  local  value 

f 

of  the  drag  coefficient.  Thus  the  quantity 


increases  from  its  initial  value  corresponding  to  (C  ),  to  the  final 

f 

value  far  downstream  corresponding  to  the  total  drag  coefficient  of  the 
body.  The  opposite  limiting  case  of  "slowly-adjusting^  turbulence  is 
characterized  by  "freezing"  the  quantity  (^f/(L)  (  ” — g—  )  YT 

00  f 

at  its  initial  value  corresponding  to  (CD  • 


A  precise  value  of  (CD  )i  is  difficult  to  determine  a  priori,  but 

a  rough  estimate  is  obtained  by  equating  (Df ^  to  the  sum.  of  the  skin 

friction  drag  of  the  body  and  the  additional  momentum  defect  produced 

by  the  pressure  rise  associated  with  the  flow  deflection  at  the  neck. 

The  analysis  in  Appendix  4  indicates  that  (C^  )±  varies  with  Reynolds 
1  * 

number  like  Re""2  when  the  free  shear  layer  is  laminar.  According  to 

Eq.  (A.  4-9)  at  Red  =  105,  (Cp  ^  ^  0.02  for  axially- symmetric  flow, 

so  that  (C*  ).  ^  0.04  at  =  6,  while  (CJ  )1  ^  0.2  at  =  20. 

1  00  f 

*  The  "exact"  expression  for  €  T  (Eq.  (13)  or  (l6))  is  utilised 
in  the  detailed  computations  of  Sections  3  and  4. 
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Tiie  actual  estimated  value  of  (C_  1 ).  is  less  significant  than  the  fact 

f 

that  it  is  at  least  one  order  of  magnitude  lower  than  the  total  body- 
drag  coeffcient  for  a  bluff  body  at  hypersonic  speeds.  According  to 


Eq,  (l8)  the  initial  value  of  — ^ — g-  Y_  (or  the  "frozen"  value  for 
"slow"  turbulence)  is  also  one  order  of  magnitude  lower  than  the  value 
corresponding  to  turbulent  flow  across  the  major  portion  of  the  outer 
wake . 


Ey  employing  Eqs.  (15)  and  (17)  we  have 


^®Y  ■ 
VA 


(  )foo  ~'j  fi  cpj 


Thus,  once  (C„  ' ).  is  known  the  initial  value  of  Y  is  determined  by 
ut  1  if 

the  value  of  B^/H  selected  to  start  the  computations.  In  other  words, 

B^/H,  (Cd  )i  ,  and  H  are  the  parameters  of  the  problem.  If  (Y^  )^«.  1 

one  expects  the  subsequent  development  of  the  wake  to  be  relatively 

insensitive  to  the  particular  value  of  (Y„  ),  selected  (Section  k). 

1 


2.3.2.  Growth  of  Inner  Weke  for  "Locaily-Similar"  Turbulence 


"S.  q» 

_ — _ i  ^  _ _ 


Although  small,  the  initial  value  of 


plays  an  important 


role  in  determining  the  initial  rate  of  growth  of  the  inner  wake  once 
the  pressure  approaches  ambient.  For  "locally  similar"  turbulence  this 
initial  growth  rate  determines  the  location  downstream  of  the  body  at 
which  the  inner  wake  begins  to  swallow  appreciable  momentum  defect  from 
the  outer  wake,  and  thereby  to  increase  - — ~ — ^ — 

CO 


in  accordance  with 
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with  Eq.  (18).  In  Figure  4,  the  momentum  defect  contained  in  the  inner 
wake  is  represented  by  the  area  under  the  curve  of  enthalpy  ratio 


hT-hi 


vs.  Y. 


(Eqs.  (12),  (l4),  and  (15)  ).  From  this  figure  one  sees  that 


'  *  *>T\ 


a  _a_ 

V  K. 


*  (20) 


where  the  right-hand  side  represents  the  swallowing  rate.  By  employing 
Eqs.  (17)  and  (20)  we  find  that 

I  dj  f )  Ka^yjoT  Tt\-v, 


and  by  using  Eqs.  (2)  and  (18)  we  have 


4 — (Si-  ±L_  Y  Ko^ttv  H  ft.  JY,  ,  (22) 

But  the  energy  balance  along  the  wake  axis  (Section  3)  leads  to  a 


relation  of  the  form 


'  ^ '  *-4 '  d  %  vnrwy 


which  together  with  Eq.  (22),  yields  a  differential  equation  for 
€.  Tj  m 

— — t -  Y-,  of  the  form 

co  f 


A. — . f  P4*  €.-r  v  rn~\ 
d^Vfcc  ‘t4/ 


r  i  /,'~ 

lw  N  \rtV-,  \  (  ^-T  v  \ 

1  ,T* 


*  A  proof  of  this  relation  is  given  in  Section  3.1. 
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This  last  relation  clearly  shows  the  exponential  character  of  the  growth 
of  the  turbulent  diffusivity,  but  the  amplification  factor  depends 
strongly  on  the  slope  of  the  inviscid  enthalpy  profile  g'(zf),  as 
expected  (Section  4). 


In  the  initial  stages  z„  4.C  1  and  g*  (z^)<<  1,  so  that  the  rate 

X  x 

of  growth  of  the  inner  wake  is  given  by  an  expression  closely  resembling 
the  usual  similarity  relation  (Section  3): 


-m.+v 

f 


_  K 

d  J  -  V"*1 


•  (23) 


By  taking  g'(zf)  =  -  2zf  (Appendix  1),  and  by  using  Eq.  (23)  and  the 
approximate  relation  (^f / f°°)  (Appendix  2),  the 

(£.  m 

differential  equation  Eq.  (22)  for  - -r —  can  be  integrated,  and 

U  CO  a 

the  result  is  as  follows,  for  1 ; 


TP--*'  3 

+  .... 


(24) 


where 


A  -  »  .  L  K<S(o) 


*m+ 


tp  +i 


and  )  denotes  a  suitable  average  value  of  (^f / ^ ) .  If  the 
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turbulent  mean  enthalpy  profile  is  chosen  to  be  a  parabola  then 
Gx  =  2/3,  G2  =  1/4,  G"(0)  =  -2  . 

Also, 

V\  j_ 

( u„-')  n/  ~  s 

over  a  wide  range  of  Mach  numbers  for  equilibrium  conditions 

-  2 

(Appendix  1).  As  an  illustration  we  take  g(z)  =  e~z  ,  so  that 

(Eq.  (A. 1-10)),  g1  =  \Fny2  and  gg  =  •5-;  then 

^  2  ^  O.36  ando(  2  ^1.35  for  C  =  1  (Eq.  (A. 1-9)). 

o 

Wiuh  K  =  0.04  (Appendix  3  )  and  =  0.9  one  obtains 

A0^  0.0033  (^/f„)3/2  (0„  )/  (85) 

0.03  (^/p„)2/3  (CDj)V3 

Fortunately  the  axial  distance  required  to  increase  the  effective 
turbulent  diffusivity  by  a  specified  factor  is  relatively  insensitive 

to  the  value  of  (c^  ). .  This  distance  varies  like  (Cn  )~  ^  in  the 

1  ,  1/4  * 

tvTo-dimensional  case,  and  like  (C_  ).  in  the  axially-symmetric 

f  -l 

case.  However,  this  distance  varies  directly  with  (^/p^)  for 

_  _i 

two-dimensional  flow,  and  like  (^/j^  for  axially-symmetric  flow. 

At  Re^  =  10  ,  (0^  )^  =  0.014  in  the  two-dimensional  case  and 
(Cn  ).  *  0.02  for  axially-symmetric  flow  (Appendix  4), 

t  _  ^ 

At  M^=  6,  (pf/fe«)  =  O.39.  According  to  Eqs.  (24)  and  (25),  the 
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m 


distance  required  to  increase  the  initial  diffusivity  /u^  d)  Y^, 
by  50  per  cent  is  about  300  diameters  in  the  two-dimensional,  case 
and  about  35  diameters  in  the  axially-symmetric  case.  At  M  =  20, 

CX> 

(o  fftoo  )'=  0.033*  Here  the  corresponding  axial  distance  is  about 
3500  diameters  ( 1 )*  for  two-dimensional,  flow  and  about  120  diameters 
for  axially-symmetric  flow.  The  marked  difference  between  two- 
dimensional  and  axially-symmetric  flow  is  caused  by  two  related 
factors:  (i)  the  hot  gas  in  the  inviscid  flow  it  concentrated  much 
closer  to  the  wake  axis  in  the  axially-symmetric  case;  (ii)  the  axially- 
symmetric  inner  wake  swallows  mass  (and  momentum  defect)  much  more 
rapidly  than  its  two-dimensional  counterpart,  again  because  of  the 
geometric  factor.  The  growth  of  ^/u^  d  for  the  axially-symmetric 
inner  wake  is  indicated  schematically  in  the  accompanying  sketch. 


*  In  other  words  this  distance  is  so  large  that  the  quantity 
(ff/f°°  )  (£rp/ua>  d)  is  virtually  constant  at  its  initial  value. 
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Mass  Concentration 


3.1.  Turbulent  Diffusion  of  Enthalpy  In  the  Inner  Wake 


In  a  steady,  two-dimensional  or  axially- symmetric  turbulent  flow 


of  "boundary  layer  type  the  energy  equation  for  the  complete  static 

7 

enthalpy  of  an  effective  binary  mixture  takes  the  following  form  : 


Here  €z  ...  and  are  the  turbulent  thermal,  momentum,  and 

mass  diffusivities,  respectively,  and  Le^  is  the  turbulent  Lewis 
number  =  The  laminar  contribution  to  the  transport  and 

dissipation  terms  is  regarded  as  negligible  in  comparison  to  the 
main  turbulent  contributions. 


The  relative  importance  of  the  dissipation  term  in  the  wake  is 
measured  by  the  ratio 


where  &u  »  uf  -  u  (0).  In  hypersonic  wakes  this  term  is  negligible. 
The  statement  that  =  1  eliminates  the  last  term  in  Eq.  (26). 

(This  statement  is  not  necessary  in  formulating  the  energy  integral 
equation.  (See  below.))  As  explained  in  Section  2  only  the  thermo¬ 
dynamic  effect  of  the  static  pressure  variation,  but  not  its  dynamic 
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effect,  is  taken  into  account  here.  With  these  approximations  Eq.  (26) 
becomes 

+  v  ■%) '  ^-(f  '-i  (27) 

12 

By  the  usual  order  of  magnitude  considerations  for  wake  flows 
v(  bh/b  y)  u  (bh/d  x)  not  too  close  to  the  neck.  Also 
u  '=  u  =  [2>  u  ,  where  0.80  =  (B  ^0.90.  Thus  the  present  problem  is 
analogous  to  a  problem  in  one -dimensional  (m=  0)  or  cylindrical 
(m  =  l)  non-steady  heat  conduction,  with  t  =  x/u.  However,  the 
situation  is  somewhat  unusual,  because  itself  depends  on  the 
amplitude  and  breadth  of  the  enthalpy  pulse  (Eqs.  (8)  and  (13)),  and 
"heat"  is  spreading  into  a  region  (y  >  yf )  containing  a  time -independent 
non-uniform  enthalpy  distribution,  but  in  which  is  effectively  zero. 
The  diffusion  of  a  material  property  is  governed  by  an  equation  similar 
to  Eq.  (27)  (Section  3-2). 


By  introducing  the  Howarth-Dorodnit syn  variable  (Eq.  (4)),  and  by 
making  use  of  the  continuity  equation,  Eq.  (27)  is  converted  to  the 
form 


(28) 
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in  Section  2  the  effective  thermal  diffusivity  £  is  defined  hy 


(f) 


, — * 

(>) 

TT 


2^ 


LL  d 

CD 


pr°)  \ 


-  (-C2L) 


€.  T 


tu  d 

CO 


Instead  of  seeking  "exact"  analytic  solutions  to  Eqs.  (27)  or 
(28)  we  integrate  Eq.  (27)  across  the  inner  wake,  thereby  converting 
it  into  a  first-order,  non-linear  ordinary  differential  equation  for 
the  enthalpy  excess,  or  momentum  defect  contained  in  the  inner  wake. 

A  second  equation  of  this  type  is  supplied,  by  satisfying  Eq.  (27 )  all 
along  the  wake  axis.**  These  two  equations  are  sufficient  to  determine 
the  behavior  of  the  two  parameters  B,  (<C)  and  (  ^ )  (Eq.  (3)),  and 
thus  all  characteristics  of  the  turbulent  diffusion  of  enthalpy. 


*  Actually  the  transformation  to  an  "equivalent  incompressible 
flow"  is  strictly  valid  only  along  the  wake  axis,  because  P„  is  not 
constant . 

**  Additional  Integrals,  or  higher  "moments"  of  Eq.  (27)  can 
be  employed  later  if  greater  accuracy  is  found  to  be  necessary. 
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Ey  multiplying  Eq.  (27)  by  ym,  integrating  across  the  inner  wake, 
and  utilizing  the  continuity  equation  and  Eq.  (4),  one  obtains 


where  u  has  been  replaced  by  u,  and  thereby  eliminated  from  Eq.  (30)* 
The  integral  on  the  left-hand  side  of  Eq.  (30)  is  proportional  to 


*  The  heat  transfer  rate  q  (yf )  at  the  interface  between  the 
inner  and  outer  wakes  (y  =  yf )  is  supposed  to  be  negligible  compared 
to  the  rate  at  which  enthalpy  is  entrained  at  this  interface.  Since 
•  /  4  /  ,  -K.  /  c)ku  \ 

=  cpL  {  d”M 

this  approximation  is  correct  if 


_  ^  «  . 


or  if 


±4-  4^  >?  1 

v£  dx 


the  momentum  defect  (Eqs.  (12)  and  (1*0 ).  By  introducing  Eq.  (3)  for 
the  turbulent  enthalpy  profile,  and  employing  Eq.  (2),  one  has 


This  last  equation  is  identical  to  the  relation  obtained  from  Figure 
4  by  intuitive  considerations  (Eq.  (20);  now  the  basic  assumptions 
underlying  this  relation  are  clear.  Integration  of  Eq.  (31)  yields 


tn+i 

Yt+ 


(32) 


taJ-l 

rM 

KcoJo  ] 

$ 

L 


"rrv-H 


By  recognizing  that  (Appendix  2) 


*m.+l 


VpJ  v- 


Tf 


a  V  ™'  -  <f  (constant)  ,  *  (A. 2-6) 
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TO4-1 

and  employing  Eq.  (A.  1-9)  for  ,  Eq.  (32)  is  rewritten  in  a 

convenient  form  defining  a  "drag  function"  F^z^): 


m-H 


ft  .  Y 

^7  ^  Tt  c 


Tn+i 


B>  Yt+ 


(33) 


+ 


£x>o  I 


(Cd4)l  tm+l 


/  hL 
I  I  T 


\ 


2  ^(^)d2  - 


"|**  ) 

ib _ iiBt) 

~rn  -vi 

SH-J 


where  zf  =c*-m  and  H  =  (hj/h^  )Q  -  1. 


Along  the  wake  axis  Eq.  (27)  or  (28)  becomes 


TV 


(34) 


( 


By  using  Eq.  (3)  for  the  turbulent  enthalpy  profile,  one  obtains 


\  G  "(°)  Cl  f  ±Z  ) 


(35) 


r4 


For  "locally  similar"  turbulence  £  is  given  by  Eq.  (l6).  Ey  utilizing 
Eqs .  (17a),  (33),  and  (A. 2-6),  one  finds  that 


( 
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£i  .  °C  ,  pc°) F'(^i) 

(3u~d  \  yC*r-=C"‘4>  (36) 


Now  (dB^/cL^  )  can  be  eliminated  in  favor  of  (dzf/dj2  )  by  employing 
Eqs.  (31)  and  (A. 2-6),  and  itself  is  expressed  in  terms  of  zf  by 
using  Eq.  (33)*  For  simplicity  we  take  P~(p/h).  Finally  one  has 


3:('2-f,H)(<iVdf)''(i</^')J7^'  ^C7cD  )i(Pf /P-f^  (37) 


Tt\-*2. 


where 


F1(zf)  is  defined  by  Eq.  (33)  5 


Tw-Vl  f 


H  °DU 

w.r  c^d 


(A. 2-8) 
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Ho) 

!■ 


t  -h  M 


(37c) 


with  (Eq.  (33)), 

P_  (.Cf)L  jr.„  M*»)  (37d) 

'  ^  /  '  *>  * .  ■"  ■*  r 


Also 


/. 


^-4l 


4 


(37e) 


(37* ) 


The  function  jL-(zf,  *0  in  Eq.  (37)  depends  on  ^  only  through  the 
slow  variation  of  H  with  decreasing  (p^/p^ )  along  the  wake  axis. 
Therefore,  (zf,  H)  is  practically  a  function  of  zf  alone,  and  Eq.  (37) 
can  he  converted  into  a  single  quadrature  for  the  growth  of  the 
turbulent  inner  wake: 


4^3 


r 

"T,  J-| 


/ 


\ 


*; . 


(38) 
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The  static  pressure  distribution  along  the  vake  axis  must  be  obtained 
from  a  theoretical  solution  for  the  inviscid  flow  including  the  wake 
shock  (Figure  l),  or  from  experiments.  Of  course  for  ^  ?  100, 

(Pf/P <»  )  *  1.  Once  zf(  ^  )  is  known  from  Eg..  (38)  the  physical 
wake  boundary  is  determined  from  Eq.  (37e  )>  and  the  amplitude  of  the 
turbulent  pulse  is  given  by  Eq.  (37<i)* 


In  the  opposite  limiting  case  of  "frozen"  diffusivity  from  Eq.  (l8a) 


(18a) 


and  by  using  Eqs.  (31),  (33 )}  and  (A. 2-6)  the  heat  conduction  equation 
(Eq.  (35))  is  converted  to  the  following  single  quadrature: 


where 


(39) 


(39a) 


*  The  exact  value  of  the  initial  diffusivity  given  by  Eq.  (l6) 
or  Eq.  (36)  is  used  in  the  calculations  discussed  in  Section  4. 
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with  F2(zf)  defined  by  Eq.  (3Tb )j.  and  (^/h^  )  given  by  Eq.  (37f)> 


In  the  region  not  too  close  to  the  neck,  but  before  appreciable 

"swallowing"  has  begun, 

'  '*■'  (  ■n'+i  ^  ✓  .  _  , 

^  A  «.  ^  <0  ,  F,(^+)=  I  ,  Fk(2+;»l  j 


KM 


\  +  H; 


C<~rrv 

-  ' 


Thus  the  solution  of  both  Eqs.  (37)  and  (39)  is  reduced  to  the  similarity 
relation 


yt4  -  (YT|). 


^  ^“rtyv* 


However,  the  outward  growth  of  the  physical  boundary  of  the  inner  wake 
is  strongly  influenced  by  the  falling  static  pressure  along  the  wake 
axis  in  this  region  (Eqs.  (A. 2-6)  and  (37e)). 

Once  the  swallowing  process  begins,  the  solutions  for  "locally 
similar"  turbulence  and  "frozen"  diffusivity  diverge,  of  course,  and 
it  is  difficult  to  give  any  precise  estimates  of  the  wake  growth.  For 
locally  similar,  axially- symmetric  turbulent  wakes,  however,  the 
swallowing  process  is  virtually  completed  100-200  body  diameters  down- 


3-11 


A  *  (°DA 


stream  of  the  neck,  so  that  F  1(z;£.}'^'  — F  2(zf)  =  1 
for  f  •’  V100  -  200.  Roughly  speaking  (Eqs.  (37),  (37a),  (3?c)  -  (3?f)) 


'-r 


) 


(h)  <™r=^  ^  '/( 


(41) 


and 


Thus  the  wake  width  at  a  given  < is  not  very  sensitive  to  H  (or  flight 
velocity).  The  growth  of  the  wake  is  given  more  accurately  by  Eq.  (37e); 
this  relation  takes  into  account  the  cooling  effect,  which  makes  the 
wake  grow  more  slowly  than  an  isothermal  wake .  Numerical  calculations 
are  required  to  illustrate  the  actual  behavior  of  the  wake.  (Section  4). 


Far  downstream  (  ^ ,>  10  )  the  wake  has  cooled  off  sufficiently 


so  that 


h  (0) 


=~1, 


h 


r~S ,  .  o(m  ^f 

=  1  ,  and  - 

+  rw 


«  1.  In  this 


'<X>  “co  "f 

regime  the  growth  of  the  wake  is  described  by  a  relation  like  Eq.  (40), 

with  C_  replaced  by  (C_  +  Cn  )  for  "locally  similar"  turbulence. 

\  0  \ 

The  enthalpy  distribution  takes  on  the  self-similar  behavior  analogous 
to  Townsend's  case: 


v  Q  / 
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3.2.  Turbulent  Mass  Diffusion  la  the  Inner  Wake 
Once  the  history  of  the  turbulent  diffuslvity ^  ( v)  along  the 
wake  axis  is  known,  turbulent  diffusion  of  a  particular  chemical 
species  in  the  inner  wake  can  also  be  described  by  an  analysis 
similar  to  that  employed  for  the  complete  static  enthalpy  (Section 
3-1).  As  an  illustrative  example,  consider  the  simplest  case  in  which 
all  of  the  diffusing  species  is  confined  within  the  inner  wake,  and 
recombination  or  other  chemical  reactions  involving  this  species  are 
so  much  slower  than  diffusion  that  they  can  be  ignored.  Such  a 
situation  might  arise  when  an  ablating  body  deposits  a  "foreign" 
substance  into  the  boundary  layer. 


f 


In  analogy  to  Eg..  (3)  for  the  enthalpy,  the  simplest  representation 
of  the  mass  concentration  is  a  two-parameter  description  of  the  form 


/ 
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where  K.  is  the  mass  fraction  of  the  diffusing  species,  and  Y_  (K  ) 
1  iD  ^ 

is  a  measure  of  the  breadth  of  the  mass  fraction  profile .  In  this 

special  case  it  is  convenient  to  take 


VT 


Y 


(*6a) 


Then, 


*Tt\ 

-*,-2(17)  pu.K,xj  p>p4  ^ 


TH-M 


i  inti 

d  Ka\  P. 

•>.  O  I  ps 


a  COTLst.  (44) 


Tfv+V 


where 


CO 


Xt'Yx  ;  r,  =  tf%  ,  Fi  *  >£ 


By  introducing  a  mass  flux  coefficient,  c^,  one  has 


ft. 


TTv-Vl 


2C-nT6f. 


^  J"*'  ,=  C'"' 

oo  a>  Y(\-*| 


—  (^COn_std_Tl.t^ 


(44a) 


A  second  relation  between  K  and  Y^  is  supplied  by  the  binary 

ip 

diffusion  equation,  which  for  le^,  -  1  takes  the  same  form  as  Eq.  (27) 
or  (28),  with  h  replaced  by  Kj.*  Along  the  wake  axis  (Eq.  (34)), 


(45) 


In  spite  of  this  fact  the  normalized  mass  fraction  and  enthalpy 
are  not  similar  during  the  swallowing  process  (Section  3)>  beca\ise  they 
do  not  satxsxy  x>xm  &£Uu0  u u uj.iu.cl jry  condi.  t i.on.3 . 


*  - 

% 
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or  (Eq.  (43a)) 


rSW 

d  § 


c 


(45a) 


For  "locally  similar"  turbulence 


•T 


pi, 


is  given  by  Eq.  (36 )j  by 


introducing  this  relation  and  Eq.  (44a)  into  Eq.  (45a)  one  obtains 


K 


Tn. 


kCCd*); 


dS  c^/<r^'0 


Yrv-H 


1  fa  \ 

fv  )  -H 


mM1*6) 

^  “S*  ) 


By  again  taking  0 — <  p/h  and  integrating,  one  finds 


l-n\ 


,  -rn+i  -m-M  / « *  ■+Tr'- 

[Zj  -  OW,  d)  £4, 

3^,n) 
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When  the  static  pressure  (p^P^  )  is  known  and  zf(%)  is  obtained  from 
the  quadrature  in  Eq.  (38),  all  quantities  in  L  (|  )  are  also  known 
(Eq.  (37)).  An  alternative  form  of  (Eq.  (47))  is  given  by 


where  the  initial  turbulent  mass  fraction  and  enthalpy  profiles  are 

2  2 

"matched"  by  taking  (Y^  ).  =  ■§•  (Y^  ).  (Appendix  3). 

XD  1  f 

Once  the  mass  fraction  is  known  the  value  of  the  number  density 
per  unit  volume  along  the  axis  is  calculated  frcm  the  relation 


The  number  density  distribution  across  the  wake  is  given  by 


where  Y„  is  obtained  from  Eq.  (48) . 
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In  the  opposite  limiting  case  of  "frozen"  diffusivity 


4 


is  given  hy  Eq.  (l8a),  with  CL  =  (CL  ).,  and  Eq.  (45a)  becomes 

vf  uf  1 

(Eq.  (44a)) 

Tr  /fln-ti)  -T^/  \ 
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so  that 


1 


In  this  case  )  is  obtained  from  the  quadrature  for  frozen 

diffusivity  (Eq.  (39) ),  and  the  number  density  per  unit  volume  is 
calculated  from  Eqs.  (49),  (50),  and  (53),  with  the  appropriate 
values  of  h/h  . 

According  to  Eq.  (49)  the  initial  decay  of  number  density 
along  the  wake  axis  is  caused  mainly  by  the  volume  expansion  of  the 
flow  as  pf  — r  p ^  .  When  (  if  -  %.  i )  x  50  -  100,  pf  =  p^  ,  and  the 
cooling  of  the  wake  flow,  or  its  increasing  density,  means  naturally 
that  the  number  density  falls  off  less  rapidly  with  downstream  distance 
than  the  mass  concentration. 


As  stated  before,  the  swallowing  process  for  the  axially- 
symmetric  turbulent  inner  wake  is  virtually  completed  when  ( t.  -  ^ ^  > 
200,  so  that  Fg  (zf)  =  1  when  (^-  ^)>200.  Also 


+  1  . 

c 


AND 


By  comparing  Eqs.  (37)  and  (47),  one  sees  that 


2  2 
Y  =  Y 

D  xf 


This  slow  decay  of  peak  concentration  along  the  axis  represents 
the  limiting  case  of  zero  or  "slow"  recombination  rate.  If  the 
diffusing  species  is  reactive,  the  equilibrium  concentration  ("very  fast" 
reaction  rates)  is  a  function  of  the  temperature  and  pressure,  and  is 
known  once  the  enthalpy  field  is  known  (Section  3*1* )•  The  actual 
distribution  of  the  diffusing  species  along  the  wake  axis  lies  between 
these  two  extremes,  and  the  methods  of  the  present  paper  can  be 
generalized  to  include  finite  recombination  rates  whenever  these  rates 
are  known. 
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4.  Distribution  of  Some  Observables  in  the  Wake :  Comparison  with 
Experiment 

4.1.  Inviscid  Enthalpy  Pro^'.l^  for  Axlally-Symmetrle  Plow  and 

Initial  Conditions  for  Inner  Wake 

Two  examples  have  "been  calculated  In  detail  for  axially-symmetric 
flow:  (l)  M  =  8.5  and  pw  =  1  ATM,  corresponding  to  the  conditions 
for  most  of  the  experimental  data  on  wake  growth  obtained  in  ballistic 
ranges;  (2)  Moa=  22  at  an  altitude  of  100,000  feet,  corre spending  to 
typical  ballistic  reentry  conditions..  The  calculations  were  divided 
into  three  different  regions,  beginning  at  the  wake  neck  and  moving 
downstream.  First,  a  very  short  region  was  assumed  to  occur  in  which 
pure  expansion  (no  mixing)  takes  place  because  of  the  rapid  drop  in 
pressure  immediately  behind  the  wake  shock.  This  effect  is  so  pre¬ 
dominant  that  it  overrides  any  mixing  that  could  occur  until  the 
pressure  is  3-^  times  ambient.  In  that  region  the  entire  flow  field 
was  expanded  simply  by  the  volume  ratio  required  to  reduce  the  neck 
pressure  to  the  initial  mixing  pressure,  which  was  assumed  to  occur 
at  p  =  hp^ .  (The  fact  that  pressure  heavily  dominates  at  first  and 
then  mixing  gradually  begins  to  take  hold  is  brought  out  clearly  by 
the  wake  growth  results.)  The  second  region  extends  from  the  station 
where  p  =  kp^to  the  point  at  which  the  pressure  becomes  ambient.  In 
this  region  both  pressure  drop  and  mixing  occur.  In  the  third  and 
final  region  the  pressure  is  ambient  everywhere  and  mixing  alone 


needs  to  be  considered. 


T’ 
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It  is  necessary  to  determine  the  neck  pressure  and  fix  its  decay 
to  ambient  conditions.  Characteristics  calculations  were  carried  out 
at  M  =  22  for  the  geometry  shown  in  the  sketch  below  using  the 
Wohlwill  computer  program.  Briefly,  the  pressure  pg  (see  sketch) 
immediately  behind  the  wake  shock,  which  depends  on  Mach  number,  is 
about  24p  for  M  =  22  (p~  ^  7*5  P  for  M  =8.5),  and  the  pressure 

(jO  &'■  c.  c&  OP 

drops  to  a  value  four  times  ambient  in  a  distance  of  about  six  body 

diameters  behind  the  wake  shock  (only  one  diameter  is  required  at 

M  =  8.5).  Since  y_  --  .18  for  these  conditions,  the  turbulent 

neck 

core  expands  by  a  factor 


to  a  value ry  =  .44  in  this  initial  region  before  mixing  begins 


(y~  =  .25  at  =  8.5).  Assuming  the  neck  position  to  be  about 
1 


CO 


1.5  "body  diameters  'behind  the  nose,  the  initial  stations  measured  from 

the  nose  where  mixing  calculations  were  begun  were  x/d  =  7*5  for 

Mq3  =  22  and  x/d  =2.5  for  M  =8.5.  The  distance  downstream  to  the 

point  where  the  pressure  returns  to  ambient  is  also  a  function  of  the 

free  stream  Mach  number.  For  M  =8.5;  (x/d)  ==  25  and  for 

P  = 

M  =  22  (x/d)  _  ~  50  . 

r’J  '  9  ~ 

At  the  initiation  of  mixing,  the  inviscid  enthalpy  profile  in  the 
outer  wake  iimst  he  specified.  Since  the  enthalpy  is  relatively  in¬ 
sensitive  to  pressure  in  the  range  1  =  p/p  =  4  the  enthalpy  profile 
corresponding  to  ambient  pressure  is  employed  in  the  calculations. 

To  find  the  inviscid  enthalpy  profile  at  p  =  p^  the  shape  of  the  how 
shock  vae  used  as  described  in  Appendix  1.  If  the  shape  is  known 
analytically  and  ^  is  approximated,  an  analytical  expression  for  the 
downstream  profile  can  be  found,  as  shewn  hy  Iykoudis111-.  Alternatively, 
if  the  entropy  distribution  behind  the  shock  is  found,  the  profile  at 

p  =  p^  can  be  determined  numerically  hy  continuity  methods,  as 

15 

described  hy  Goulard  .  In  the  present  paper  the  continuity  method 
was  used  at  M  «  8.5  and  both  methods  were  employed  for  M  =  22. 

The  resulting  profiles  when  transformed  to  the  Kowarth- Dor odnit syn 
variable  hy  Eq.  (l)  are  shown  in  Figure,:  5  and  6. 

At  the  axis  the  slope  of  the  profile  must  he  zero,  as  shown  in 
Appendix  1.  Since  the  numerical  and  analytical  shock  shapes  used 
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for  =  22  all  agree  fairly  veil  in  the  tody  nose  region  (Appendix  1, 
Figure  19),  the  enthalpy  profiles  also  agree  quite  veil  near  the  axis. 
Away  from  the  axis  the  parabolic  shock  g-'ves  a  higher  enthalpy  because 
the  slope  of  this  shock  surface  is  somewhat  steeper.  For  computational 
purposes  the  enthalpy  profiles  were  approximated  by  combining  a 
Gaussian  curve  near  the  axis  with  an  exponential  curve  in  the  outer 
portion,  in  such  a  manner  that  the  total  drag  coefficient  was  unity, 
and  the  slopes  of  the  two  curves  agreed  at  the  matching  point.  It 
should  be  emphasized  that  this  procedure  was  employed  for  convenience 
on3.v,  and  exact  values  from  the  profiles  could  be  used  just  as  well. 

It  was  found  from  the  computer  program  described  previously  that  the 
wake  shock  raises  the  value  of  the  enthalpy  along  the  stagnation  stream 
line  by  5  per  cent,  Combining  these  features  gives  the  following 
equations  for  the  enthalpy,  which  are  also  shown  in  Figures  5  and  6. 
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Based  on  the  estimates  described  in  Appendix  4,  the  values  of  the 
initial  drag  contained  in  the  turbulent  core  were  taken  as  follows: 
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(y 

8.5 

.02? 

22.0 

.0077 

The  effect  of  different  values  of  (C^  )i  on  the  resultant  vake  growth 
is  considered  later  for  one  of  the  examples.  The  experimental 
constant,  K,  used  in  the  diffusivity  was  taken  as  0.04  (Appendix  3)* 
Also,  the  value  of  J5>  (Eq.*  12)  was  taken  to  he  0.9* 

4.2.  Turbulent  Diffusivity 

The  growth  of  the  turbulent  diffusivity  for  "locally  similar" 

turbulence  has  been  discussed  in  a  general  way  in  Section  2;  we  now 

examine  the  numerical  results  for  «„=  8.5.  The  exact  expression 

for  <=■„  Y_  m  can  be  obtained  from  Eqs.  (13)  and  (33)  and  is  given 
if 
by 


This  quantity  compared  to  its  initial  value  is  shown  in  Figure  7  as 
a  function  of  downstream  distance  from  the  initial  mixing  station. 
The  rapid  initial  growth  in  Y^  ra  is  caused  by  the  pressure 
decrease  and  accompanying  rapid  increase  in  y^.  The  region 
x/d1^  10  -  20  illustrates  the  eurly  stages  of  mixing  before  the 
diffusivity  increases  rapidly  again.  The  explosive  growth  of 
diffusivity  begins  at  x/d  ^  40  for  the  condition  shown.  For  two- 


4-6 


dimensional  flow  and  the  same  external  conditions  this  growth  would 
occur  much  later.  According  to  Figure  7  the  process  of  swallowing 
enthalpy  or  momentum  defect  originally  contained  in  the  outer  wake  is 
virtually  completed  in  about  300  body  diameters.  Far  downstream  the 
diffusivity  approaches  the  asymptotic  value  given  by  the  relation 


eT  v  ^ 

P  T' 
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<6x. 
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Cr 


4.3.  Growth  of  the  Turbulent  Inner  Wake 

The  calculated  growth  of  the  turbulent  inner  wake  at  a  flight 
velocity  of  9>500  ft ./sec.  is  shown  in  Figure  8  for  both  "locally 
similar"  turbulence  and  "frozen"  diffusivity.  Up  to  x/d  =  15-20  the 
two  limiting  cases  are  practically  identical,  because  the  marked 
increase  in  turbulent  diffusivity  in  the  "locally  similar"  case  has 
not  yet  begun  (Figure  7),  and  the  pressure  drop  plays  an  important 
role.  Downstream  of  this  region  the  "locally  similar"  wake  grows 
rapidly  to  a  width  of  about  3  body  diameters  at  x/d  =  100,  5  body 
diameters  at  x/d  =  300,  and  8  body  diameters  at  x/d  =  1,000.  Soon 
after  the  swallowing  process  is  completed  (x/d^300)  the  wake  width 
approaches  the  familiar  l/3-power  growth  curve  corresponding  to  the 
total  drag  of  the  body.  Wake  widths  calculated  on  the  assumption  of 
"frozen"  diffusivity  are  smaller  by  a  factor  of  at  least  two  over  the 
range  200<.x/d  ^,4,000. 


These  calculations  are  compared  with  experimental  shadowgraph  data 
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in  Figure  9.  Most  of  the  data  are  supplied  "by  Slattery1^  and  Clay 
(9,000  ft ./sec.)  at  the  Lincoln  Laboratory  (M.I.T.).  Some  additional 
points  obtained  in  the  region  x/dv^300  were  taken  from  Feldman  (AVCO) 
(15,000  ft. /sec.)  and  Dana^and  Short  (CONVAIR)  (7>400  ft. /sec.). 

The  experimental  scatter  is  too  large  to  detect  any  details  of  a 
rapid  growth  region  in  the  early  stages  (x/d<^50).  The  most  important 
result  is  that  the  predicted  growth  of  the  inner  wake  for  "locally 
similar"  turbulence  is  in  good  agreement  with  the  shadowgraph  measure¬ 
ments  in  the  range  200  \  x/d<(4,000,  while  the  wake  widths  calculated 
for  "frozen"  diffusivity  are  too  low  by  a  factor  of  about  two.  No 
adjustment  of  the  diffusivity  constant  K  was  necessary  to  fit  the 
data. 


The  theory  (Section  3)  indicates  a  weak  dependence  of  wake  width 
on  flight  velocity,  and  agrees  qualitatively  with  the  fact  that 
Feldman's  data  obtained  at  15,000  ft. /sec.  lies  in  general  above  the 
computed  curve  (Figure  9)>  vhile  the  data  of  Dana  and  Short  at 
7,400  ft. /sec.  lie  below  the  curve. 

At  a  flight  velocity  of  22., 000  ft  ./sec.  the  calculations  show 
that  the  inner  wake  grows  more  rapidly  in  the  region  x/d<300  than 
at  9,500  ft. /sec.  (Figures  10  and  11),  because  of  the  combined  effect 
of  a  larger  pressure  drop  and  a  higher  effective  turbulent  diffusivity. 
However,  the  wake  grows  more  slowly  than  the  1/3 -power  during  some 


region  x/d.,^300  -  400,  because  the  core  is  still  much  hotter  than 
ambient  (Eq.  (37®)  )•  The  enthalpies  along  the  axis  and  the  front 
ere  shown  in  Figure  12.  At  ^  =  1,000  the  axis  enthalpy  is  still 
over  six  times  ambient.  In  fact,  the  calculated  curve  of  wake 
growth  has  not  yet  taken  on  the  l/3 -power  behavior  at  x/d  =  4,000. 
Of  course  very  far  downstream  the  wake  boundaries  for  all  flight 
velocities  must  eventually  coincide  (for  the  examples  shown  this 
occurs  at  -^>10')  if  the  total  body  drag  coefficient  is  the  same. 


4.4.  Effect  of  (C  ). 

_ Jf 

In  order  to  determine  how  accurately  the  value  of  (C^  )i  must 
be  known,  calculations  were  made  at  =  22  for  two  values  of 
initial  drag  differing  by  a  factor  of  two,  namely  .0077  and  .0155. 
The  enthalpy  increment  on  the  axis  of  the  inner  wake  was  taken  to 
be  the  same  for  both  values  of  (C_  ). .  According  to  Eqs.  (A. 2-8) 

X  1 

and  (37d),  if 

(B-l/h).0077  =  (Bl/H'.0155 


then  the  quantity 

fo/H)= 
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is  also  the  same.  In  taking  =  const,  we  must  have  the  ratio 


m+1 


.Z  \ 

Tcrrr 


=  constant;  thus  if  (C_  ).  is  doubled  z_ 

x  1  ti 


must  increase 
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by  VT.  Then  the  initial  fronts  are  not  quite  at  the  same  position  and 
h  / 0  a 

neither  — — -  or  -r -  have  the  same  values.  This  discussion 

hco  heo 

shows  that  the  value  of  z„  must  be  chosen  to  be  consistent  with 

ri 

(CD  )±  and  (B1/H)i  . 

The  growth  of  the  turbulent  inner  wake  is  shown  in  Figure  13  for 
initial  values  of  the  drag  coefficient  equal  to  .0077  and  .0155  with 
and  without  pressure  effects.  The  axial  coordinate  is  measured  from 
the  initial  mixing  station.  Differences  in  downstream  location  for 
a  given  wake  width  are  large  in  the  initial  stages  (x/d  <,100 ),  but 
soon  become  unimportant  once  the  swallowing  process  is  in  full  swing. 
For  these  two  cases,  the  values  of  the  integrand  of  Eq.  (38),  J'  (zf,  H), 
were  essentially  identical  after  x/d  =  125 .  The  differences  in  this 
function  occur  only  for  small  zf  and  upon  integration  eventually 
become  a  small  contribution. 

Here  the  effect  of  the  pressure  drop  in  shaping  the  early  portions 
of  the  front  is  very  evident,  resulting  in  a  parabolic  behavior  as 
opposed  to  the  slow  initial  growth  predicted  by  mixing  alone. 

4.5.  Mass  Diffusion  Calculations 

Calculations  for  turbulent  mass  diffusion  at  M  =22  are  shown 

<20 

in  Figure  l4.  As  described  in  Section  3-2,  this  case  represents  the 
limiting  behavior  in  which  some  species  is  contained  entirely  within 
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the  inner  wake,  and  no  recombination  occurs,  so  that  the  total  number 
of  particles  remains  the  same.  The  behavior  with  "locally  similar" 
turbulence  and  "frozen"  diffusivity  is  shown  as  computed  from  Eqs.  (48), 
(49),  and  (53  )•  A  rapid  decrease  in  concentration  occurs  first,  because 
of  the  pressure  drop,  but  by  x/d  =  50-0  the  peak  concentration  along  the 
axis  is  falling  slowly  because  of  the  slow  growth  of  the  wake  width 
and  the  cooling  effect  (Figures  10  and  11).  In  these  calculations, 
the  pressure  drop  accounts  for  a  decrease  by  about  a  factor  of  five 
while  the  mixing  accounts  for  about  two  orders  of  magnitude  up  to 
x/d  =  1,000.  The  values  of  peak  concentration  for  "frozen"  diffusivity 
are  almost  an  order  of  magnitude  above  the  values  for  "locally  similar" 
turbulence . 


4.6,  Electron  Distribution 

Electron  concentrations  in  the  wake  can  be  found  from  the  enthalpy 
and  pressure  distributions.  Representative  results  are  shown  in  Figures 
15,  1 6,  and  17.  In  these  illustrative  examples  it  is  assumed  that  the 
flow  is  in  complete  thermodynamic  equilibrium  at  least  up  to  the  neck 
location.  The  lower  curve  in  Figure  15  shows  the  axial  distribution 
of  electrons  for  M  =  22  at  100,000  ft.,  assuming  that  the  electrons 

CO 

also  follow  the  equilibrium  flow  conditions  in  the  wake  given  by  the 
"locally  similar"  turbulence  calculations.  Here  the  initial  concen¬ 
tration  is  equal  to  the  value  at  the  neck,  or  in  other  words,  it  is 
given  by  the  conditions  — =  .6  - and  p  =  24  p  . 
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The  pressure  drop  to  ambient  has  a  strong  influence  for  a  short 
distance,  after  which  -ft  —  determines  the  distribution.  The 

h  CO 

distribution  crosses  the  equivalent  C-band  (56 00  mc/sec;  ng  '= 

11.  3.*  . 

3.1  x  10  /cm  )  electron  concentration  at  x/d  =  100,  L-band 

(1300  mc/sec;  ng  =  1.7  x  10^/cm^)  at  £  200,  and  the  lowest 

frequency,  UHF  (400  mc/sec;  n  =  1.6  x  10^/ cm^),  at  x/d  =  325* 

The  curve  in  Figure  15  labelled  "pure  diffusion"  corresponds 
to  the  case  of  frozen  flow  in  the  wake  so  far  as  the  electrons  are 
concerned.  The  initial  conditions  are  the  same  as  for  the  equi¬ 
librium  case,  but  the  presence  of  electrons  in  the  outer  wake  is 
ignored.  Thus  the  normalized  distribution  of  number  density  is 
identical  with  the  diffusion  results  shown  in  Figure  l4.  Again 
the  pressure  drop  in  the  initial  region  accounts  for  a  considerable 

portion  of  the  decrease  in  n  .  In  fact  the  differences  between  the 

e 

two  limiting  cases  of  frozen  and  equilibrium  flow  are  not  large  for 
x/d <100.  Downstream  of  x/d  =  150  the  electron  concentration  is 

*  In  rationalized  units,  the  electrons  /cm0  are  given  by 
me  2  f2 

n  =  - ^ - CJ  =  - n —  where  cd  and  f  are  radar  frequencies 

6  10° 

in  radians  and  cycles  per  second  respectively,  is  the  dielectric 
constant  of  free  space,  and  e  and  m  are  the  eleBtronic  charge  and 
mass  respectively.  We  consider  here  the  lengths  to  be  roughly  the 
distance  between  the  object  and  the  position  where  the  electron 
density  falls  below  the  value  given  by  this  equation  for  any  radar 
frequency. 


4-12 


orders  of  magnitude  larger  in  the  case  of  pure  diffusion,  as  expected. 
As  discussed  in  Section  3*2,  these  two  cases  should  represent  the 
limits  of  behavior,  and  the  actual  situation  lies  between  these  two 
limits . 

The  third  curve  in  Figure  15  shows  the  axial  electron  distribu- 

2 

tion  for  the  laminar  equilibrium  condition  as  calculated  by  Feldman 
l4 

;  and  lykoudis.  The  strong  pressure  effect  initially  is  similar,  but 

j 

|  once  the  pressure  is  ambient  (  ^  /> 100)  the  slow  time  scale  compared 

|  with  turbulent  diffusion  is  evident. 

i 

I 

In  Figure  l6  the  calculations  have  been  scaled  to  an  altitude 
of  200,000  ft.  for  illustrative  purposes.  Here  the  equilibrium 
conditions  would  indicate  no  C-band  wake.  The  L-band  wake  length 
corresponds  to  a  value  of  =  50,  while  the  trail  length  for  UHF 
corresponds  to  g  =  150.  However,  the  pure  diffusion  wake  would 
have  a  UHF  length  of  x/d  ^  250. 

Finally,  the  radial  distribution  of  electrons  at  several  axial 
stations  is  shown  in  Figure  17  for  M  =22  and  h  =  100,000  ft., 

CO  CO  1  7 


*  The  laminar  calculations  do  not  account  for  the  initial 
,  heatfed  inner  core,  hence  the  difference  in  the  concentrations  for 

I  j  f  <100. 


under  the  assumption  of  equilibrium  and  "locally  similar"  turbulent 
flow  conditions.  Near  the  neck  the  electrons  appear  restricted  to 
a  fairly  narrow  region,  but  they  diffuse  outward  rapidly  downstream 
of  the  neck.  The  total  width  of  the  plasma  cylinder  corresponding 
to  the  plasma  frequency  at  L-band  remains  virtually  constant  in  the 
region  30<’x/d<^110  at  about  3*5  body  diameters.  The  width  corre¬ 
sponding  to  UHF  remains  practically  constant  at  about  5*5  body 
diameters  in  the  region  30<sx/d<v250. 
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5.  Conclusions  and  Future  Work 

1.  At  velocities  of  the  order  of  9>000  -  10,000  ft ./sec.,  the 
predicted  growth  of  the  inner  vake  for  "locdly  similar"  turbulence 
is  in  good  agreement  vith  shadowgraph  measurements  of  wakes  behind 
spheres  obtained  in  ballistic  ranges.  Wake  widths  calculated  for 
the  opposite  limiting  case  of  a  "frozen"  diffusivity  determined 
only  by  the  initial  drag  of  the  inner  vake  are  too  small  by  a  factor 
of  at  least  tvo  over  the  range  200 <,  x/d  < k,  000  covered  by  the  experi¬ 
mental  data.  We  conclude  that  the  adjustment  time  for  the  scale  and 
intensity  of  turbulence  must  indeed  be  short  compared  to  the  time 
scale  for  the  changes  in  turbulent  diffusivity  along  the  wake  axis. 
The  adjustment  is  probably  made  easier  by  the  fact  that  the  process 
of  swallowing  enthalpy  or  momentum  defect  originally  contained  in  the 
outer  wake  is  almost  completed  in  about  300  body  diameters.  The 
constant  K  appearing  in  the  expression  for  turbulent  diffusivity  in 
the  present  paper  corresponds  directly  to  Townsend's  "universal" 
Reynolds  number  for  low-speed  turbulent  wakes.  Tentatively  one 
concludes  that  the  turbulence  mechanism  in  the  wake  with  respect  to 
a  fixed  observer  is  probably  quite  similar  to  the  low  speed  case,  in 
spite  of  the  large  mean  temperature  gradients.  Experimental  studies 
of  the  structure  of  turbulence  in  the  wake  are  required  to  settle 
this  point. 

Measurements  of  the  growth  of  the  turbulent  inner  vake  behind 
blunt-nosed  bodies  of  various  nose  biuntness  ratios  at  velocities 
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of  the  order  of  20,000  ft ./sec.  vould  help  to  clarify  the  influence  of 
the  nose  drag,  the  outer  vake  enthalpy  distribution,  and  the  hot  core 
in  the  inner  vake. 

2.  Once  the  swalloving  process  is  in  full  swing  the  influence  of 

the  initial  drag  coefficient  (CL  ).  is  small,  because  the  drag  absorbed 

r  1 

from  the  outer  inviscid  vake  is  so  much  larger  than  the  initial  drag. 

However,  the  "incubation  length"  required  for  the  inner  wake  to  begin 

swallowing  momentum  defect  at  an  appreciable*  rate  is  proportional  to 

-1/3  -±/k 

(C_  ).  in  the  two-dimensional  case,  and  to  (CL  ).  '  in  the 

Df  1  uf  1 

axially-symroetric  case.  A  brief  analysis  shows  that  (C^  )i  varies  with 

_i  f 

Reynolds  number  like  (Re  )"^  when  the  free  shear  layer  is  laminar,  and 

a  rough  estimate  establishes  the  value  of  (C^  within  a  factor  of  ( 

about  two.  Careful  theoretical  and  experimental  studies  of  the  re- 

compression  just  behind  the  body  are  required  if  the  initial  conditions 

for  the  turbulent  inner  wake  are  to  be  specified  more  precisely. 

Theoretical  investigations  of  the  stability  of  laminar  flow  in  the 
wake,  and  extension  of  hot-wire  studies  of  laminar- turbulent  transition 
in  the  axially-symmetric  wake  vould  be  very  helpful  in  determining  the 
lower  Reynolds  number  limit  for  the  appearance  of  turbulence. 

3.  Because  the  fluid  in  the  free  shear  layer  above  the  zero  stream¬ 
line  is  decelerated  and  compressed  at  the  "neck",  the  gas  in  the  inner 


vake  is  much  hotter  initially  than  the  gas  in  the  outer  vake  at  the 
turbulent  front.  In  one  typical  example  (M  =22,  h  =  100,000  ft.) 
the  enthalpy  along  the  vake  axis  requires  a  distance  of  about  150 
body  diameters  to  cool  dovn  to  a  value  equal  to  the  initial  enthalpy 
in  the  outer  vake  at  the  front.  Even  at  x/d  =  1000  the  enthalpy  on 
the  vake  axis  is  still  about  10  times  ambient  for  a  blunt  body  such 
as  a  sphere.  The  axial  variation  of  the  enthalpy  at  the  front  (y  =  yf ) 
is  simply  a  one-to-one  mapping  of  the  curve  of  vake  grovth  y  ( <.? )  on 
the  original  enthalpy  distribution  in  the  outer  vake . 

4-.  As  an  example  of  turbulent  mass  diffusion  the  decay  of  peak, 
electron  density  along  the  axis  of  the  inner  vake  is  calcv bated  for 
M  =  22,  h  =  100,000  ft.,  on  the  assumption  that  the  initial  electron 

JO 

density  corresponds  to  thermodynamic  equilibrium,  and  that  the  electron- 
ion  recombination  rate  in  the  vake  is  negligibly  small.  The  peak 
electron  density  decays  from  a  value  of  5  x  lO1^/ cm^  to  a  value  corre¬ 
sponding  to  the  plasma  frequency  at  L-band  (f  =  1300  me/sec; 
ng  =  1.7  x  10  /cnr)  in  a  dovnstream  distance  of  about  1000  body 
diameters.  This  length  should  be  compared  vith  the  corresponding 
distance  of  about  200  body  diameters  for  complete  thermodynamic  equi¬ 
librium  (very  fast  recombination  rates).  The  actual  situation  is 
expected  to  lie  somevhere  betveen  these  tvo  limits.  In  any  event  the 
predicted  radar  trail  length  for  a  blunt  body  is  fairly  long,  even  vith 
turbulent  diffusion. 


One  Interesting  result  for  equilibrium  flov  at  M  =  22  and 
h  =  100,000  ft.  is  that  the  width  of  the  plasma  cylinder  correspond- 
to  the  plasma  frequency  at  L-band  remains  virtually  constant  at 
yifcw.vt  3  .5  body  diameters  in  the  region  30-C  x/d-CHO,  while  the  width 
corresponding  to  UHP  (400  mc/secj  1.6  x  10  ,/cm'3 )  remains  practically 
constant  at  about  5*5  body  diameters  in  the  region  3CK.x/d<v250. 

Similar  calculations  can  be  carried  out  for  equilibrium  gas 
radiation  intensity  or  radiation  from  an  ablating  species  in  the 
wake. 


5.  The  results  obtained  in  the  present  paper  are  sufficiently 

encouraging  so  that  one  can  begin  to  think  of  including  the  effects 

of  chemical  and  electron-ion  recombination  processes.  Clearly  these 

rate  processes  and  turbulent  diffusion  interact  when  streamlines 

enter  the  spreading  inner  wake,  because  of  the  effect  of  diffusion 

on  the  local  mass  fractions  of  the  reacting  species.  Hopefully, 

18 

the  methods  utilized  by  S.  C.  Lin  and  J.  D.  Teare  (and  others) 
for  the  outer  wake  can  be  combined  with  the  present  analysis  to  give 
a  more  complete  picture  of  the  wake  at  hypersonic  speeds. 
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APPENDIX  I 

ENTHALPY  DISTRIBUTION  IN  THE  OUTER  WAKE 

x 

We  consider  the  inviscid  flow  region  for  ^  >  25-50  where  the  static 

pressure  is  virtually  equal  to  the  ambient  pressure.  The  enthalpy 

profile  is  independent  of  axial  distance  and  is  specified  by  the  bow 

and  wake  shock  shapes.  The  procedure  used  here  to  calculate  the 

profile  is  to  determine  the  profile  that  would  exist  for  a  semi-infinite 

hemisphere -cylinder,  and  then  modify  this  profile  slightly  to  correct 

2 

for  the  wake  shock.  Feldman  has  shown  that  this  correction  is  small 
and  for  many  applications  could  be  neglected. 

If  the  bow  shock  shape  is  known  analytically  and  ^  is  approximated, 

an  analytical  expression  for  the  downstream  profile  at  ambient  pressure 

14 

can  be  found,  as  shown  by  Lykoudis ,  Alternatively,  given  a  numerical 
shock  shape,  the  mass  flux  and  entropy  distributions  just  downstream 
of  the  shock  can  be  found.  Assuming  an  isentropic  expansion  to  ambient 
pressure  some  distance  downstream  of  the  nose,  the  flow  variables  can 
be  found  corresponding  to  the  ambient  pressure  and  the  entropy  values 
downstream  of  the  shock.  The  radial  position  or  distribution  of  the 
streamlines  is  determined  by  a  continuity  balance  for  the  mass  flux 
entering  the  bow  shock  along  any  given  streamline.  Goulard1'’  first 
applied  this  method  as  a  means  of  comparing  various  flow  field 
techniques . 
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Thus  the  problem  is  basically  that  of  determining  the  correct  shock 
2 

shape.  Feldman  used  the  method  of  characteristics  and  computed  down- 

X 

stream  to  ^  =  25  along  the  hemisphere -cylinder.  These  values  were 

19 

subsequently  used  by  himself  and  the  Goulards  to  solve  the  laminar 
heat  conduction  problem  in  an  equilibrium  wake  flow.  However,  as 
Goulard^  pointed  out,  the  characteristics  method  is  subject  to 
accumulative  errors  when  it  is  extended  to  the  downstream  distances 
needed.  The  present  authors  have  had  the  same  experience  with 
computer  calculations  run  in  conjunction  with  this  investigation. 
Considerable  care  in  the  programming  procedure  is  necessary. 

Blast  wave  theory  indicates  that  a  parabolic  shock  should  exist. 
lk 

Iykoudis  has  used  this  fact  with  the  decided  advantage  that  a  closed 
solution  was  obtained  even  to  the  laminar  conduction  problem.  There 
are  two  fe attires  which  modify  this  assumption:  (l)  Because  of  the 

entropy  layer,  a  parabolic  shock  can  not  exist  for  a  hemisphere- 

•* 

cylinder  and  the  exponent  must  be  slightly  less  than  .5-  This  result 
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has  also  been  found  by  Van  Hise  ,  who  used  the  sonic  cone  method  and 
obtained  computer  solutions  to  various  body  shapes.  He  found  the 
shock  shapes  for  all  the  various  bodies  could  be  correlated  very  well 
by  the  equation 


(A.l-l) 

(2)  Real  gas  effects  (the  Van  Hise  computer  calculations  were  for 
constant  )  can  be  important .  The  density  no  longer  has  the  given 

*This  fact  was  recently  pointed  out  to  us  by  Dr  Milton  Van  Dyke 
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maximum  ratio  -  across  the  shock,  and  the  shock  shape  is  modified 

$  “  * 

accordingly. 

Regarding  the  correction  for  the  secondary  shock,  it  has  been  found 
that  the  body  streamline  enthalpy  when  expanded  to  p^_  is  raised  by 
about  5$  over  that  on  the  hemisphere-cylinder  surfa-e.  The  amount  in 
axially- symmetric  flow  can  be  somewhat  different  from  that  for  two- 
dimensional  flow,  since  the  free  expansion  streamline  behind  the  body 
in  the  former  case  tends  to  curve  inward  as  the  axis  is  approached, 
resulting  in  a  different  secondary  shock  strength. 

For  the  specific  examples  used  in  the  present  work,  the  experi¬ 
mental  shock  shapes  shown  in  Figure  18  were  employed  to  get  the 
entropy  and  mass  distribution  for  M  =8.5.  The  entropy  distribution 

CO 
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downstream  of  this  shock  was  found  using  the  Wohlwill  computer 

shock  routine,  and  then  the  continuity  balance  described  before  was 

used  to  get  the  profile  at  p  =  p  .  At  Mw  =  22  the  Wohlwill1^ 

computer  program  was  used  to  get  the  subsonic  and  transonic  shock 

shapes.  The  result  at  =  22  is  shown  in  Figure  19  compared  with 

20  l4 

the  shock  shapes  given  by  the  Van  Hise  and  Iykoudis  equations. 

For  both  examples,  the  final  enthalpy  profiles  when  transformed  by 

the  Howa.'fbh-Dorodnitsyn  variable  are  shown  in  Figures  5  and  6.  | 


/ 


The  laminar  profile  exhibits  a  flat  region  near  the  axis  followed 
by  a  rather  rapid  drop  in  enthalpy.  Since  this  shape  contributes  to 


A-1.4 


the  initially  slow  end  then  explosive  growth  of  the  turbulent  mixing, 

one  should  examine  theoretically  whether  the  initial  slope  at  the 

axis  is  finite  or  indeed  zero.  This  point  can  be  investigated  by 

21 

using  the  procedure  given  by  Wohlwill  for  calculating  the  gradients 
normal  to  a  streamline  crossing  a  shock.  We  first  calculate  the 
entropy  gradient. 


"T  d  3  __  -y  5  dep 
3  dip  3  dep  dip 


(A. 1.2) 


where  \  is  the  stream  function,  Tg  is  the  temperature  behind  the 
shock  and  cj)  is  the  shock  inclination  arigle  measured  from  the  axis. 
Since 


T5dS  =  A--#- 


(A-1.3) 


and  from  the  Rankine  -  Hugoniot  relations 


k.  =  K.^  +  i  <4.  [>  -(-jr9  ] 

t>  *  p.,,  sL-twK'--|ri)+  p» 


it  follows  that 


(A.  1.4) 


(A. 1.5 ) 


(A. 1.5a) 
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We  want 
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For  two  -  dime  ns  i  onal  flow, 


»C  5in4co5<j>(|--j^9  ii>. 
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dR  = 


pudr 


and  for  axially- symmetric  flow, 

^>»w-00R-=1R  =  fu-^dr 

where  the  notation  is  shown  in  the  sketch. 
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Now  is  a  constant  along  any  streamline.  Since  the  pressure 
constant,  TdS  =  dh  and  we  have  finally  for  two-dimensional  flow 

Td5  =  dk  =  Tfpu-juj  cxx$Q-$rf 
dv'dr  Ts(f„“te;R\s 

and  axially- symmetric  flow 

Tdl=_dk=  Tfipujr  u.t  (  >--^r 
dT  dr  Ts  (p~  V)  Ra 

Since  cj?  =  90°  for  the  stagnation  streamline,  the  normal  entropy 


(A. 1.5b) 


(A. 1.6a) 


(A. 1.6b) 


is 


(A. l.?a) 


(A. 1.7b) 

and 


enthalpy  gradients  are  always  zero  for  two-dimensional  flow  regardless 
of  the  body  shape.  For  axially-symmetric  flow,  the  gradient  is 
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proportional  to  tbs  "body  radius  r,  and  is  zero  only  if  the  streamline 
returns  to  the  axis,  which  it  does  in  the  present  case.  Thus  the 
enthalpy  gradient  at  the  axis  is  zero. 


Finally,  with  the  enthalpy  profile  of  the  form  given  hy  Eq.  (2) 

fe  - ')- =  fe  1 1  (°0'u )  -  )  »> 

the  value  of can  he  found  from  the  drag  coefficient  .  When 
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the  drag  coefficient  is  given  hy  the  expression 
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APPENDIX  2 

RELATION  BETWEEN  INVISCID  AND  TURBULENT 
MAPS  FLUX  QUANTITIES  YL  “+1  AND  (pf  /  f>w  )  Y^  m+1 


The  rate  of  entrainment  of  mass  flux  by  the  inner  wake  is 
given  by  the  expression 
dm.  dy 

sr-  =  ?f  uf ( ar  -|5f)  (a-2-1) 

where  is  the  inclination  of  the  streamline  just  entering  the  inner 

wake.  In  the  present  analysis  the  influence  of  the  turbulent  inner 

wake  on  the  outer  wake  is  regarded  as  a  second-orde'v  effect  (Section  2). 

In  other  words,  uf,  andjb^  are  supposed  to  have  the  same  values 

along  the  curve  y  =  yf(x)  as  they  .would  have  if  the  outer  wake  flow 

field  extended  smoothly  to  the  axis.  To  this  approximation  Eq.  (A.2-1) 

states  that  the  rate  of  increase  of  mass  flux  m^.  contained  within  the 

bounding  curve  y  =  yf(x)  is  also  the  same;  therefore,  the  inviscid  and 

turbulent  mass  flux  quantitities  at  any  station  differ  by  a  constant, 

* 

at  most. 


This  constant  is  readily  determined  by  evaluating  the  two  mass 
flux  quantities  at  x  =  x^  where  Y^  m+1  and  (^f/ )  Y^  m+1  are  both 


a*. 


*Near  the  neck  (dy^/dx)  and  (2>f  may  be  comparable,  and  some  error 
could  be  introduced  bynegiectingthe  Interaction  effect.  However,  the 
wake  flow  is  "pressure-controlled"  in  the  first  5-10  body  diameters 
downstream  of  the  neck  in  any  case.  Beyond  x/d  =  10,  «  (dy^/dx) . 
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<<.  1.  When  zf  «.l,  1  +  H  =  (|vu/go  ,  and  according  to 

Eq.  (37d), 

c^ru\frY'+H) 


(A. 2-2) 


or 


~(YLf)  (A, 2-3) 

By  inverting  Eq.  (4),  taking  ^f/p=  ,  and  utilizing  Eq.  (3),  one 

finds  on  the  other  hand  that 

By  employing  Eq.  (19),  we  have 
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and  "by  equating  the  two  relations  Eqs.  (A. 2-3)  and  (A. 2-5)  for  y^,  one 


obtains 
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where  (Eqs.  (A. 1-9)  and  (A.2-6)) 
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APPENDIX  3 

SOME  REMARKS  ON  TOWNSEND'S  EXPERIMENTAL  RESULTS  FOR  LOW- 
SPEED  WAKES:  RELATION  BETWEEN  OUR  K  AND  TOWNSEND'S  R^ 


According  to  Townsend's  experimental  results  X 


near  the  axis  of  a  symmetrical  low-speed  wake  behind  a  long  cylindrical 


rod  normal  to  a  uniform  air  stream,  where  £  ^  is  independent  of  y. 
Once  the  wake  flow  becomes  "self -preserving"  the  mean  velocity  distri¬ 


bution  near  the  wake  axis  is  described  by  the  parabola 


where  a  is  a  constant.  Thus 


(A.3-1) 


(A. 3-2) 


where  Au  =  uf  -  u(O).  But  T=  -  p  U.V  ,  and  the  question  arises  as 
to  whether  Eq.  (A. 3-2)  for  f  is  consistent  with  the  behavior  of  the 


correlation  u'v''  across  the  wake . 


Townsend  found  that  the  turbulent  velocity  fluctuations  require 
a  somewhat  longer  distance  behind  the  rod  to  reach  equilibrium,  or 
"similarity",  than  the  mean  axial  velocity.  Nevertheless  this  simi¬ 
larity  is  well  established  at  x/d  =  500  (Figures  7*1  and  J.h,  p.  135 
and  pp.  142-14-3,  Reference  6).  In  Section  2.1  it  was  pointed  out  that 
the  characteristic  length  in  our  case  is  not  the  body  diameter  but  the 
inner  wake  momentum  thickness,  0.  Therefore  on  the  basis  of  Townsend's 
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results  and  our  estimates  of  ©i/d,  one  should  expect  complete  flow 
similarity  in  the  inner  wake  to  he  established  in  5  -  20  body 
diameters.  Once  this  similarity  exists 


i  i 

cu  v 


id/1.-) 


and  by  symmetry 


(A. 3-3) 


(A.3-3a) 


near  the  wake  axis  (Figure  7«10,  page  152,  Reference  6).  The  quantity 

r\ 

u,<=  is  virtually  independent  of  y  near  the  wake  axis  and  its  magni- 

*  2 

tude  decays  with  downstream  distance  exactly  like  ( A^u)  (Figure  7.4, 

— j?  p 

page  142,  Reference  6).  In  other  words  u’  —  (Au)  .  By  combining 
this  statement  with  Eq.  (A.  3-3)  and  (A.  3-3a),  we  find  that 

r  s-pu.V-  - p (A  (/H  Aj i)  (A. 3-4) 


Clearly  Eq.  (A. 3-2)  and  (A. 3-4)  for  T  are  compatible  if  <=.T  -'•Au.'Ljf, 

AilU'U.l. 

or  if  the  quantity  — — — *-*-  is  a  universal  constant  "Reynolds  number" 
for  the  symmetrical  wake  flow.  In  this  particular  case,  the  validity 
of  the  semi-empirical  relation~=  p£-T(^u/c)a|J  is  assured  by  symmetry  and 
similarity;  in  fact  one  could  dispense  with  it  entirely  and  work  directly 
with  the  Reynolds  stress.  Evidently  these  same  arguments  are  applicable 
to  axially-symmetric  "similar  wakes",  and  by  Reynolds'  analogy,  to  the 
turbulent  transfer  of  mass  and  heat. 


The  only  remaining  question  is  the  magnitude  of  the  Reynolds  number 
^  ,  which  must  be  determined  by  rough  physical  arguments  about  the 
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characteristic  size  and  shape  of  the  large  eddies  ,  or  by  experiment  . 

The  procedure  adopted  here  is  to  make  certain  that  the  relation  between 
the  turbulent  diffusivity  and  the  drag  coefficient  of  the  wake  (Eq.  (l8a)) 
reduces  to  Townsend's  relation  for  isothermal  flow.  In  this  way  our  K 


is  related  directly  to  Townsend's 

t, 


(A. 3-5) 


T 


where  X  is  a  measure  of  the  wake  width.  In  Townsend's  case,  the 
o 

mean  axial  velocity  distribution  not  too  near  the  edge  of  the  wake 
is  described  very  closely  by  the  Gaussian  distribution 


O.  oc 

,J^oa  U(o) 


exp. 


Thus  (Eq.  (l4)) 
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where  F^  =  ^/Z  and  F2  =  ■§•.  By  using  Eq.  (A. 3-5)  and  (A.3-7) 
we  have 

’a  C 
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(A. 3-8) 


TT.  +  1 


Now  for  isothermal  wakes  of  Townsend's  type  our  Eq.  (l8a)  for  . 


(or 
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)  reduces  to  the  form 


u._-d 


Yrf)  = 
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(A. 3-9) 
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By  equating  Eq.  (A.  3-8)  and  (A.  3-9)  for  <==•  one  finds  that 
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(A. 3-10) 


For  a  two-dimensional  wake  with  a  parabolic  mean  enthalpy  profile, 

=  2/3,  and 

K  =  0.53  (Brjv)"1  (A.3-10a) 


or,  taking  Townsend's  value  of  R^  =  12.5; 
K  =  0.043 


(A.3-10b) 


In  the  case  of  an  axially- symmetric  wake,  the  value  of  K  depends 

to  some  extent  on  the  "matching"  between  the  velocity  profile  selected 

for  the  inner  wake  in  our  case  and  the  Gaussian  velocity  profile  for 

a  low-speed  wake  spreading  into  an  irrotational  flow.  By  comparison 

6  9 

with  the  experimental  results  of  Townsend  ,  and  Hall  and  Hislop,  the 

"best  fit"  for  a  parabolic  velocity  profile  is  obtained  by  taking 

-8(1  /d).  At  this  point  (u^  -  u)  ^  0.05  Au  for  the  actual 
Tf  o 

profile  (Figure  7.4,  Reference  6),  and  the  turbulent  intensities 
2  2  2 

u'  ,  v'  ,  and  w'  are  less  than  10  per  cent  of  their  peak  values.* 


*  The  problem  of  selecting  the  "best  fit"  is  somewhat  similar  to 
the  problem  of  choosing  the  "thickness"  of  a  boundary  layer  along  a 
solid  surface.  One  wants  to  include  the  bulk  of  the  mean  shear  flow, 
but  not  too  much  of  the  "external"  inviscid  flow.  Fortunately  in  the 
axially- symmetric  case  the  wake  width  at  a  given  downstream  location 
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With  this  choice  of  ,  and  with  Gg  =  l/k ,  and  R^  =  14.5 


K  =  0.035  =  0.04  (A.3-10c) 

This  value  of  K  is  utilized  throughout  the  numerical  computations 
(Section  4). 


*  &  depends  on  K 


1/3 


(Section  3-1); 


or  on 


However,  the 


downstream  location  for  a  given  wake  width  varies  inversely  with  K.  It 
should  he  noted  that  the  intermittency  factor  ^  is  about  0.4  at  the 
point  (Figure  7*5>  Reference  6),  and  a  sensitive  hot-wire 


anemometer  measuring  turbulent  fluctuations  might  find  a  "wake  width" 
as  much  as  50  per  cent  greater  than  the  wake  width  for  the  mean  flow 
quantities . 
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APPENDIX  4 

MOMENTUM  DEFECT  IN  THE  INNER  WAKE  AND  ORDINARY  DRAG; 

ESTIMATE  OF  INITIAL  DRAG  COEFFICIENT,  (C  ). 

Ut 

According  to  Eq.  (l8)  the  effective  turbulent  diffusivity  is 
determined  by  the  momentum  defect  D^  in  the  inner  wake  measured 
with  respect  to  the  local  flow  in  the  outer  wake  at  the  turbulent 
front.  But  the  usual  aerodynamic  drag  D^  is  measured  with  respect 
to  the  free-stream  velocity,  u  .  These  two  quantities  are  con- 

CO 

nected  by  the  relation 

=  Zip)'™  f  Pt  = 

0  (A. 4-1) 


or, 

s  £>oc  "  (A. 4-2) 

where  m^  is  the  mass  flux  in  the  turbulent  inner  wake .  Even  in 
the  initial  stages  near  the  neck,  where  D^  is  very  nearly  constant, 

D  increases  in  the  downstream  direction  because  the  inner  wake  is 
capturing  additional  mass  flux,  and  the  term  m^,  (u^  -  u^,)  takes 
account  of  the  differences  in  momentum  flux  as  measured  in  the  two 
different  coordinate  systems. 
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Once  the  process  of  swallowing  momentum  defect  begins  and 
D  both  vary  along  the  wake  axis.  Of  course  the  total  momentum 
defect  is  constant  when  (  P  -  £«,)  80  *s  'btle 

difference  between  the  turbulent  and  inviscid  drags  of  the  inner 
wake  for  the  same  mass  flux.  This  relation  can  be  derived  as 
follows : 

According  to  Eq.  (A. 4-2)  the  difference  ADf  can  be  written 
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(A. 4-4) 


(A. 4-5) 


The  first  Integral  in  Eq.  (A. 4-5)  represents  the  drag  Dq  of  the  original 
flow  in  the  outer  wake,  and  is  a  constant;  therefore  the  quantity  in 
|  j  ,  which  is  Identical  with  AD^,  is  also  a  constant.  The  constant 
is  determined  by  evaluating  A Df  when  yf  1.  In  that  case  uf. 


(  / 


and 
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AD,  -  ait) 


-m 


pr  U“T  CU7  U"T)  ^  ; 


which  is  exactly  the  initial  momentum  defect  of  the  inner  wake  with 
respect  to  the  local  flow  in  the  outer  wake  at  the  front.  The  fact 
that  ADf  is  constant  can  also  he  verified  by  utilizing  Eq.  (20).; 
the  rate  of  increase  of  the  first  term  in  Eq.  (A. 4-3),  which 
determines  the  turbulent  diffusivity,  is  exactly  balanced  by  the 
rate  of  increase  of  the  second  term. 

A  precise  value  of  (Df )^  or  (C^  i)  is  difficult  to  determine  at 
present,  but  a  rough  estimate  is  obtained  by  equating  (Df )  to  the 
sum  of  the  skin  friction  drag  on  the  body  and  the  additional  mo¬ 
mentum  defect  produced  by  the  pressure  rise  at  the  neck  (Figure  l). 
According  to  the  momentum  equation,  this  increase  in  momentum  defect 
is  given  by 
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Pu/ar 


-m  "rTVH 


(A. 4-6) 


where  Ap  is  the  pressure  rise  at  the  neck.  Viscous  stresses  are 
neglected  because  the  compression  occurs  over  a  short  distance  . 
Let  the  subscripts  1  and  2  denote  conditions  just  upstream  and 
downstream  of  the  neck,  respectively.  Then  Ap  =  p>,  f —J^  — 
and  (Eq.  (A. 4-6)) 


For  two-dimensional  flow  the  quantities  pg  and  p^/pco  can  he 
evaluated  approximately.  C.  F.  Dewey,  Jr.* **  points  out  that  for  an 
isentropic  compression  process  p^  is  equal  to  the  pressure  on  the 
body  surface  at  the  point  where  the  streamline  is  parallel  to  the 
free  stream  direction.  By  employing  a  Prandtl-Meyer  expansion 
from  the  sonic  point  on  a  blunt-nosed  body  this  pressure  is  evalu¬ 
ated,  and  by  using  Chapman’s  criterion  for  Pg/p^  Dewey  finds  that 
p-j/p^  ='  1.3  x  10“2  (  ^  M  2).  Also  P2/P1  3  because  M1  ~  3, 

virtually  independent  of  M 

CO 

The  remaining  unknown  is  the  neck  width  (y^/d)^  This  width 
is  determined  by  the  boundary  layer  thickness  at  the  separation 
point  on  the  body  surface,  plus  the  subsequent  mass  flux  entrainment 
in  the  free  shear  layer.  A  rough  estimate  of  (y^/d^  for  the  case 
of  a  circular  cylinder  is  obtained  by  considering  this  whole  process 
to  be  equivalent  to  the  development  of  the  laminar  boundary  layer 
over  a  blunt-nosed  slab  of  equivalent  length.  Since  the  free  shear 
layer  makes  an  angle  of  about  17°  with  the  free  stream  direction 
the  neck  location  corresponds  approximately  to  =2.3.  Calcu¬ 
lations  similar  to  those  made  in  Reference  22  for  a  hemisphere- 
cylinder  yield  (cj*/d)i  '/ie^  ^18,  or(^/d>^ 


*  Private  communication . 

**  The  boundary  layer  in  the  two-dimensional  case  is  thicker  than 
it  is  on  the  hemisphere-cylinder  by  the  factor  \[2. 
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This  rough  estimate  is  in  fair  agreement  with  measurements  of  the  neck 
width  in  the  wake  behind  a  circular  cylinder  made  by  C.  P.  Dewey,  Jr., 
at  M  =  5«8* 

GO 


In  the  axially- symmetric  case  the  inviscid  flow  is  more  complicated 
near  the  neck.  For  the  present,  we  utilize  the  estimates  for  p^/p^ 
and  Pg/P^  obtained  from  two-dimensional  flow  considerations.  If  the 
ratio  of  mass  flux  in  the  inner  layer  at  the  neck  to  the  mass  flux  at 
the  boundary  layer  separation  point  is  the  same  as  in  the  two-dimensional 
case,  then  mass  continuity  dictates  that  the  quantity  (y^/d)^  must  be 
equal  to  the  two-dimensional  value  of  (y.j>/d).,  divided  by  the  factor  tf~2. 

Collecting  all  of  these  estimates  we  find  that 


v  _  "Vz 
—  2 
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(A. 4-8) 


Now  CD  =  Cf  +ACd  ,  where  Cf  is  the  skin-friction  drag  coefficient 

f±  f  ^  (/ .3) 

of  the  body.  Roughly  C„  7  — 2 — —  } 

so  that 


(A. 4-9) 


Since  the  probable  minimum  Reynolds  number  for  turbulent  flow  in  the 

3  4 

wake  is  of  the  order  of  5  x  10  ,  the  maximum  relevant  value  of  (C  ) 

D_r.  i 
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is  about  0.02  in  the  two-dimensional  case,  and  about  0.028  in  the 
axially-symmetric  case.  This  condition  corresponds  to  a  flight 
altitude  of  about  l8o, 000  feet  for  a  body  one  foot  in  diameter  at 
satellite  velocity.  At  an  altitude  of  100,000  feet,  for  example, 
Re,  =  2,5  x  10^  for  the  same  body,  so  that  (CT.  ).  ^ 0.0033  in  the 

i. 

two-dimensional  case  and  (CL.  ).  =  0.005  in  the  axially-symmetric 

f 

case.  In  the  numerical  examples  of  Section  4  the  value  of 

0.022  for  (CD  )i  at  M^=  8.5  corresponds  to  a  Reynolds  number  near 

the  lower  limit,  while  the  value  0.0077  for  M^=  22  corresponds  to 

Re ,  =  8  x  10  . 
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Figure  1.  Wake  Behind  Blunt  Body  at  Hypersonic  Speeds 


(a)  UPSTREAM 


WAKE  AXIS 


Figure  Z,  Typical  Inner  Wake  Velocity  and  Enthalpy  Profiles  just 
Upstream  and  Downstream  of  Neck 


Figure  6.  Outer  Inviscid  Enthalpy  Profile  (M^  =  22  and  100,000  ft  Altitude) 
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igure  8.  Turbulent  Wake  Width  (M^  =  8.5  and  p  =  1  ATM) 


Figure  10.  Turbulent  Wake  Width  (M  -  22  and  100,000  ft  Altitude) 


Fig.,  m  11  Cor-tpari  .on  oi  Tut  bulent  Wake  Width 


Figure  IE.  Effect  of  Initial  Turbulent  Drag  Coefficient  on  Wake  Growth 


(.gu-re  13*  Enthalpy  Di^tr ibutior..  Along  the  Axial  Streamline 


Figure  14.  T  irbulent  Diffusion  in  the  Wake 
(M^  =  22  and  100,000  ft  Altitude) 


M.  *  22  ,  ht.  ■=  100,000  FT, 

"locally  similar"  turbulence 


'C  BANtV' 


-LAMINAR  EQUILIBRIUM, 


LYKOU 


/ASSUMING  ZERO  f 
RECOMBINATION,  PURI 
^DIFFUSION 1 | — 


L  BAND* 


/ASSUMING  ELECTRONS 
"  FOLLOW  EQUILIBRIUM  , 
-  FLOW  CONDITIONS  — f 


mm 

■I 


Figure  15.  Electron  Densities  Along  the  Axial  Streamline 
{M^  =  22  and  100,000  ft  Altitude) 
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Figure  16.  Electron  Densities  Along  the  Axial  Streamline 
(M^  =  22,  100,000  and  200,000  ft  Altitudes) 


Figure  17.  Radial  Electron  Density  Distribution 
(h'l  =  22  and  100,000  ft  Altitude) 
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CHARACTERISTICS,  STL 
COMPUTER  RESULTS 


Figure  19.  Bow  Shock  Shape*  for  a  Hemisphere  at  M  =  22 
100,000  ft  Altitude  00 


